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To 

Multiply  by 

mils 

millimeters 

0.0254 

inches 

centimeters 

2.54 

feet 

meters 

0.3048 

ml  lea 

kilometers 

1.6093 

square  inches 

square  centimeters 

6.4516 

square  feet 

square  meters 

0.0929 

square  miles 

square  meters 

2,589,998.0 

cubic  inches 

cubic  centimeters 

16.38706 

cubic  feet 

cubic  meters 

0.0283 

cubic  yards 

cubic  meters 

0.764555 

gallons  (U.S.) 

liters 

3.785 

gallons  (Imperial) 

liters 

4.542 

ounces 

grams 

28.349 

pounds 

kilograms 

0.454 

tons  (short) 

kilograms 

907.185  i 

tons  (long) 

kilograms 

1,016.047 

pounds  per  foot 

newtons  per  meter 

14.59390 

pounds  per  square  inch 

newtons  per  square 
centimeter 

0.6894757 

pounds  per  cubic  inch 

kilograms  per  cubic 
centimeter 

27,679.90 

pounds  per  square  foot 

newtons  per  square 
meter 

47.88026 

pounds  per  cubic  foot 

kilograms  per  cubic 
meter 

16.0185 

Inches  per  second 

centimeters  per  second 

2.54 

inch-pounds 

meter-newtons 

0.1129848 

inch-kips 

meter-kilonewtons 

0.0001129848 

Fahrenheit  degrees 

Celsius  degrees  or 
Kelvins® 

kilotons 

tera joules 

4.183 

axo  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings* 
use  C - (5/9) (F- 32).  To  obtain  Kelvin  (K)  readings*  use  K ■ (5/9) 

(F— 32)  + 273.15. 
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1. 


INTRODUCTION 


A methodology  for  systematically  evaluating  the  uncertainties  associated  with 
the  prediction  of  vulnerability  and  hardness  of  ground  facilities  has  been  used  for  the 
past  ten  years  to  develop  confidence  statements,  from  the  defender's  point  of  view,  for 
weapon  systems  survivability  prediction.  In  this  application,  the  defender  theoretically 
lias  perfect  intelligence  of  site  conditions  and  structural  configurations  ani  materials, 
and  can  establish,  analytically  or  experimentally,  the  weapon  effect  environments  and 
structural  fragilities  to  within  acceptable  uncertainty  bounds.  Parametric  sensitivity 
analyses  are  conducted  to  evaluate  the  impact  of  various  uncertainties  on  the  system  sur- 
vivability. This  information  proves  to  be  Invaluable  for  conducting  design  trade-offs, 
and  for  weapon  effects  research  planning. 

The  present  study  represents  the  first  attempt  to  apply  such  a methodology  to 
analyzing  the  uncertainties  in  target  vulnerability  and  hardness  predictions.  Target 
vulnerability  and  liardness  and  associated  uncertainties,  given  good  intelligence,  were 
evaluated  for  two  representative  targets:  (1)  an  above-surface  industrial  building;  and 
(2)  a shul low-bur  led  command  and  control  bunker. 

The  procedure  for  evaluating  target  hardness  uncertainty  was  started  with  the 
formulation  of  mathematical  models  of  the  weapon  effects  and  the  target  response  and 
failure.  Structural  responses  models  for  both  targets  were  developed  based  on  classical 
structural  dynamics  considerations.  However,  for  the  buried  banker,  the  structural  model 
is  more  uncertain  because  of  the  inadequate  data  base.  Uncertainties  on  model  parameters 
were  characterized  and  propagated  through  the  analysis  using  a Tay 1 or-ser les  expansion 
method,  to  evaluate  the  best  estimate  hardness  and  the  hardness  uncertainty  as  a function 
of  weapon  yield.  The  results  were  evaluated  and  compared  with  the  Physical  Vuluet ab  1 1 1 1 y 
System. 

Uncertainties  caused  by  potential  systematic  errors  of  analysis,  l.e.,  bias 
between  the  mathematical  model  and  the  real  world,  were  distinguished  from  nonsystemut ic , 
random  variations  of  physical  properties.  When  propaguted  through  the  weapon  effect 
scaling  laws,  nonsystemat Ic  variations  lead  to  variations  in  estimated  damage  distances. 
These  nonsvstemat lc  damage  distance  variations  are  treated  adequately  by  the  Physical  Vul 
nerabillty  System,  and  for  most  cases  do  not  cause  significant  errors  in  the  estimated 
damage  probabilities.  However,  systematic  errors  in  target  hardness  estimates  cause 
systematic  errors  In  estimated  damage  probabilities  which  are  not  treated  adequately  by 
the  Physical  Vulnerability  System  at  present.  The  emphasis  of  this  studv  it:  focused  upon 
the  problems  of  quantifying  the  effects  of  the  systematic  errors  of  analysis. 
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Proa  the  result*.  It  wee  concluded  that  for  above  surface  targets  such  as 
Industrial  buildings,  the  current  procedure  for  determining  the  best  eetiaate  of  physical 
vulnerability  is  adequate.  Uncertainty  in  the  target  physical  vulnerability  can  be  quan- 
tified baaed  on  the  state-of-the-art  knowledge  of  low  level  blast  environments,  material 
properties,  and  structural  response  modeling.  Scale  model  tests  of  industrial  targets 
would  not  be  worthwhile  unless  th«  structural  models  were  accurate  representations  of  a 
high-value  class  of  uniformly  constructed  structures. 

For  the  shallow-buried  targets  su'h  as  command  bunkers,  procedures  for  physical 
vulnerability  prediction  to  achieve  the  desired  confidence  levels  remain  to  be  validated. 
Of  significant  uncertainty  is  the  structural  response  model  and  failure  criteria. 

Testing  of  model  structures  to  failure  with  the  objective  of  developing  more  accurate 
prediction  techniques  is  recommended.  Additional  sensitivity  analyses,  and  research  cost- 
benefit  trade-off  studies  should  be  conducted  to  provide  a basis  for  setting  priority 
for  reaaarch  projects  pertinent  to  bunker  vulnerability. 
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2.0  STUDY  APPROACH 


Target  hardness  predictions  are  generally  derived  analytically,  and  the 
uncertainties  in  these  predictions  are  dependent  upon  the  inputs  to  the  analysis,  and 
also  upon  the  analytical  model  Itself.  In  the  present  study,  mathematical  models  for 
nuclear  weapon  environ*  ats  and  target  atructural  loads  and  response  were  developed 
as  a first  step.  The  sources  of  uncertainties  in  the  model  parameters  were  then 
identified,  and  nominal  values  and  uncertainty  bounds  for  these  parameters  were 
determined.  The  analytical  models  were  used  to  determine  both  target  hardness  predictions, 
and  the  sensitivity  of  these  predictions  to  perturbations  within  the  uncertainty  bounds 
of  the  parameters  input  to  the  models.  Next,  the  uncertainties  of  the  parameters  were 
propagated  through  the  structural  response  model  to  determine  the  uncertainty  bounds  on 
the  target  hardness  estimates.  The  final  and  most  difficult  step  in  the  quantification 
of  the  target  hardness  uncertainty  was  to  estimate  the  accuracy  of  the  mathematical  model 
itself,  end  incorporate  this  information  into  the  target  hardness  estimate. 

Section  2.1  discusses  the  distinction  between  random  variations  and 
uncertainties,  and  how  these  quantities  are  defined.  The  methodology  for  calculating 
hardness  uncertainties,  and  the  Impact  of  these  uncertainties  on  targeting  is  given  in 
Section  2.2. 

2.1  STATISTICAL  CONSIDERATIONS 

If  a number  of  cylinders,  cast  from  the  same  batch  of  concrete  and  cured  in  a 
carefully  controlled  manner,  were  subjected  to  an  unconfined  compression  test  in  the 
same  loading  device  at  essentially  the  same  time,  identical  resistance  would  not  be 
exhibited  by  each  cylinder.  Individual  differences  in  the  cylinders,  much  as  air  voids 
in  the  concrete,  which  occur  randomly  and  can  never  be  completely  eliminated,  would  cause 
cylinder-to-cylinder  strength  variations.  Thus,  the  strength  of  cylinders  cast  from  a 
single  batch  of  concrete  is  a random  variable.  Random  variations  are  the  actual  physical 
variations  that  occur  within  a group,  such  as  cylinders  cast  from  the  same  batch  of 
concrete,  or  from  the  same  concrete  mix  specifications. 

When  a specific  cylinder  is  considered,  the  unconfined  compressive  strength  is 
not  a random  variable,  since  it  has  one  specific,  true  value.  This  true  value  c«n  be 
estimated,  for  example,  based  on  a knowledge  of  the  properties  of  the  concrete  batch 
from  which  it  was  cast.  However,  such  an  estimate  will  vary  from  the  true  value  of 
strength  of  this  particular  cylinder  because  of  random  variation  among  cylinders  from  the 
batch.  The  true  value  of  the  unconflned  compressive  strength  of  a particular  cylinder 
can  be  determined  by  subjecting  the  cylinder  to  a compression  test. 
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Even  after  the  test  has  been  performed,  there  remains  some  uncertainty  in  the 
result.  This  determination  of  the  true  value  of  the  compressive  strength  of  the  cylinder 
is  only  as  good  as  the  accuracy  of  the  test  will  permit.  Errors  in  the  test  measurements 
will  cause  an  error  in  the  determination  of  the  true  value  of  the  cylinder  strength. 

In  many  cases,  this  uncertainty  can  have  more  impact  on  the  outcome  of  a study  than  the 
random  variations.  Effects  of  raudom  variations  tend  to  be  self  cancelling,  while 
uncertainty  often  adds  a significant  error  to  the  entire  study  result. 

Uncertainty  refers  to  potential  differences  between  the  real  world  and  the 
perceived  world  of  analytical  models  and  of  model  parameter  values.  There  will  be  uncer- 
tainty in  the  measure  of  central  tendency  of  a random  variable  (for  example,  the  mean  or 
median  value),  in  the  measure  of  its  spread  (for  example,  the  standard  deviation),  and 
in  che  distribution  function  selected  for  the  variable.  The  uncertainty  in  any  of  these 
statistical  properties  is  the  confidence,  or  the  expectation,  that  the  true  parameter 
value  is  less  than  a specified  value,  defined  over  the  range  of  possible  values.  (A  50% 
confidence  value  Implies  that  the  actual  parameter  value  is  as  likely  to  be  higher  as  it 
is  to  be  lower) . Uncertainty  in  the  measure  of  central  tendency  of  a targeting  analysis 
variable,  or  in  the  analytical  model  chosen  could  introduce  significant  uncertainty  in 
the  targeting  analysis  results.  However,  uncertainties  in  the  variance  or  in  the 
distribution  function  of  an  analysis  variable  are  generally  not  significant. 

The  bounding  values  of  a variable,  or  the  limiting  responses  associated  with 
the  extremes  of  modal  assumptions,  are  the  usual  bases  for  estimating  uncertainties. 

There  is  necessarily  a degree  of  subjectivity  on  the  confidence  level  ascribed  to  bounds. 
Reasonable  engineering  bounds  are  sometimes  taken  to  represent  95%  confidence  bounds. 
Thus,  if  the  uncertainty  in  the  variable  under  consideration  is  normally  distributed, 
these  bounds  will  correspond  to  a region  within  plus-and-minus  two  standard  deviations 
from  the  mean. 

In  targeting  analysis  both  random  variation  and  uncertainty  effects  should 
be  properly  treated.  The  response  of  a specific  concrete  structure,  subjected  to  the 
effects  as  ^elated  with  a particular  nuclear  device  detonated  at  a specific  location  and 
time,  is  deterministic.  There  will  be  uncertainty,  but  not  random  variation,  in  the 
response.  This  is  true  even  though  some  parameters  of  the  analysis  model  might  be 
described  in  terms  of  both  a random  variation  and  an  uncertainty  distribution,  as  ior 
example,  when  the  concrete  mix  specifications  are  known,  but  evlinder  test  data  for  the 
actual  pour  are  not  available.  The  random  variance  associated  with  the  mix  strength 
actually  represents  an  uncertainty  on  the  cylinder  strength  of  the  concrete  in  the 
specific  structure  of  interest.  The  effect  of  this  uncertainty  is  combined  with  that  of 
other  sources  of  uncertainty  to  determine  the  total  uncertainty  in  the  target  hardness. 
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Where  there  are  many  similar  structures  constructed  using  the  saute  mix 
specifications,  the  random  variation  in  mix  strength  will  have  a smaller  effect  on  the 
uncertainty  in  the  median  structure  response  if  all  structures  are  subjected  to  similar 
nuclear  loading  conditions.  Random  variation  in  concrete  strength  causes  some  structures 
to  be  stronger,  and  others  to  be  weaker.  This  variation  tends  to  average  out  when  the 
median  response  of  a large  number  of  "tr'ictures  is  considered.  Thus,  when  analyzing 
many  structures,  random  variation  and  uncertainty  cannot  be  combined.  Since  the  random 
variations  tend  to  be  dominated  by  the  weapon  miss  distance  (CEP)  distribution,  their 
effects  are  generally  not  critical. 

An  appreciation  of  the  distinction  between  random  variation  and  UP'  ertainty 
leads  to  better  understanding  of  the  research  options  available,  in  addition  to  aiding 
the  understanding  of  complex  targeting  situations.  For  example,  the  ratio  of  the 
standard  deviations  of  the  uncertainty  and  the  random  variations  will  influence  the 
number  of  tests  required  to  reduce  the  uncertainty  to  an  acceptable  level. 


In  the  present  study,  random  variations  in  the  analysis  parameters  are  assumed 
to  cause  only  random  variations  in  the  system  response.  However,  it  is  possible  that  the 
effects  of  random  parameter  variations  around  the  median  conditions  will  not  average  out, 
even  after  a large  number  of  similar  structures  are  considered.  If  this  were  the  case, 
a systematic  bias  would  be  introduced  in  the  response.  For  example,  if  a group  of 
buried  structures  have  a very  stiff  median  foundation  modulus,  those  sites  with  higher 
than  median  stiffness  might  exhibit  nearly  the  same  response  as  the  median  site,  i.e.  , 
the  response  associated  with  a rigid  founiation.  However,  those  sites  with  much  lower 
stiffness  could  have  a markedly  different  response  caused  by  a rigid  body  mode.  The 
interaction  between  random  variations  and  uncertainties  can  possibly  be  removed  by 
categorizing  sites,  such  that  the  response  phenomenology  is  the  same  for  all  sites  in  a 
category,  and  by  performing  separate  targeting  analyses  for  each  category.  However, 
this  interaction  is  treated  in  more  formal  statistical  survivability  modeling,  such  as 
in  the  FAST  code  (Reference  1). 


Another  approximation  used  in  this  study  is  that  the  uncertainty  distributions 
are  assumed  to  be  lognormal..  This  assumption  in  commonly  used  in  weapon  effects  analysis 
because  ic  conveniently  approximates  much  of  the  data.  If  there  are  a number  of  uncertain 
parameters  in  the  system,  and  these  parameters  are  multiplicative  factors  affecting  the 
system  response,  then  the  Central  L.mit  Theorem  indicates  that  these  uncertainties  will 
be  lognormally  distributed  (Reference  2). 
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Consistent  with  the  assumption  of  a lognormal  distribution,  best  estimate 
values  of  the  parameters  are  taken  as  the  median  value  of  the  distribution.  For  a 
continuous  variable,  the  median  is  defined  as  the  value  of  the  variable  at  vhich  SOX 
of  the  area  under  the  distribution  function  lies  above  the  median,  and  SOX  below.  For 
a lognormal  distribution,  the  median  corresponds  to  the  mean  of  the  logarithm  of  the 
variable.  Uncertainties  are  defined  in  terms  of  2o  K-factors  which  represent  approximately 
9SX  bounds.  This  means  there  is  approximately  a 95X  expectation  that  the  true  value  of 
the  parameter  will  be  within  a factor  of  K and  1/K  of  the  median  value  of  the  parameter. 

The  K-f actor  is  related  to  the  standard  deviation  of  the  distribution  of  the  logarithm 
of  the  variable  by 

2a  - In  K (2-1) 


2.2  EVALUATION  OF  UNCERTAINTIES 

Two  primary  factors  determine  the  importance  of  the  uncertainties  in  various 
parameters  affecting  target  response:  (1)  the  sensitivity  of  the  target  hardness,  p,  to 
uncertainties  in  the  input  parameters,  x^;  and  (2)  the  uncertainty  in  the  median  values 

of  the  parameters.  If  the  target  hardness  is  written  as 

In  p • f(ln  x^)  (2-2) 


then  from  Reference  2 (Equation  4.39),  the  standard  deviation  in  the  uncertainty  of 
In  p,  Op^,  which  is  caused  by  the  uncertainty  in  the  particular  parameter,  x^ , can  be 

computed  from  the  approximate  relationship 


°Pj 


iQn.&L. 

3 (In  Xj) 


(2-3) 


where  the  partial  derivative  gives  the  sensitivity  of  the  target  hardness  to  variations 
in  xj  and  oj  is  the  standard  deviation  in  the  uncertainty  of  In  x j . The  approximation 

arises  from  neglecting  terms  higher  than  first  order  in  the  Taylor-series  expansion  of 
the  function,  f,  used  in  calculating  the  standard  deviation.  Then,  from  Equation  (2-1) 


In  K 


Pj 


P>, 

3 (In  Xj) 


m 


(2-4) 


where  K . and  K.  are  the  uncertainty  factorc  corresponding  to  o , and  o , respectively. 
PJ  J PJ  j 

Continuing  with  the  approach  of  Reference  2 (Equation  4.44),  an  expression  for  the 

total  target  hardness  uncertainty  factor,  K , can  be  derived  as 

P 
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(2-5) 


where  the  p^  are  the  correlation  coefficients  between  the  target  hardness  uncertainty 


factors. 


Equations  (2-4)  and  (2-5)  mnke  possible  the  evaluation  of  the  relative  Impact 
of  the  parameter  uncertainties  on  the  target  hardness  uncertainty.  In  addition,  the 
potential  benefits  in  reduced  target  hardness  uncertainty,  resulting  from  reducing  the 
uncertainties  in  the  key  parameters  can  be  readily  evaluated  with  these  expressions. 
Whether  these  uncertainties  can  be  reduced  by  test  must  be  determined  by  a separate 
evaluation  of  the  potential  errors  associated  with  the  available  test  techniques.  If 
no  suitable  techniques  exist,  then  the  requirements  for  simulation  capability  can  be 
stated  in  terms  of  the  noeded  uncertainty  reduction. 


To  determine  the  need  for  uncertainty  reduction,  the  target  hardness 
uncertainties  are  propegated  through  the  scaling  laws  for  weapon  effects  to  determine 
damage-distance  relationships,  and  the  corresponding  uncertainties,  as  a function  of 
yield.  Then,  the  Impact  of  these  uncertainties  on  the  yield  and  accuracy  requirements 


of  the  attacker  can  be  evaluated. 


For  this  study,  each  target  was  considered  to  be  the  object  of  a single  attack 
aimed  directly  at  that  target.  Each  target  was  one  of  many  constructed  from  the  same 
designs  and  sited  in  the  same,  generic  geology.  Therefore,  random  variations  in  as-built 
material  properties  are  about  15  to  30%,  and  construction  tolerances  are  usually  small. 
For  a given  generic  geology,  weapon  effect  variations  are  also  small,  even  considering 
the  effect  of  weather  conditions  (which  are  partially  known),  ard  total  random  variation 
on  target  damage  distances  will  usually  have  a standard  deviation  (o^)  less  that  20% 

of  the  median  damage  distance  (R) , Based  on  the  detailed  analysis  of  Reference  3, 


(Appendix  B),  the  probability  of  damage,  P , is  very  nearly  the  same  as  the  cookie  cutter 

K. 


equation: 


P »1-P 
K cc 


(2-6) 


if  Pg<0.9  and  o^<0.2R.  The  cookie  cutter  approximation  is  that  there  is  zero 


probability  of  failure  at  ranges  greater  than  R and  unity  at  ranges  clcser-in  than  R. 
Then,  convoluting  this  damage-distance  model  with  the  weapon  miss  distance  Rayleigh 
distribution  gives 


Jf 
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exp  (-0.5^/o^) 
0,849  CEP 


(2-7) 


A property  of  the  Rayleigh  distribution  is  that  there  is  a 50%  probability  of  the  weapon 
landing  within  a distance  from  the  target  equal  to  the  CEP  (circular  probable  error). 

The  random  physical  variations  do  not  significantly  affect  the  probability  of 
damage  computation  in  the  above  case.  Their  effect  ip  much  smal’er  than  that  caused  by 
random  variation  of  weapon  miss  distance  (the  CEP  distribution).  On  the  other  hand,  a 
systematic  uncertainty  in  u damage  distance  estimute  does  affect  the  confidence  iu  the 
calculated  probability  of  damage. 

The  attacker's  confidence  (expressed  as  a function  of  damage  distance)  is  the 
expectation  that  uncertainties  in  the  targeting  analysis  would  not  result  in  a damage 
distance  greater  than  u given  value.  Confidence  estimates  are  thus  based  on  the 
systematic  uncertainty  in  damage  distance,  assuming  a lognormal  distribution  of  the 
uncertainty.  The  kill  probability  as  u function  of  confidence  level  may  be  determined 
from  Equation  (2-6),  using  this  damage  distance  versus  confidence  level  distribution. 
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3*0  INDUSTRIAL  TARGET 

This  section  describes  the  analysis  of  a typical  heavy  steel  frame  industrial 
building,  as  shown  in  Figure  3-1,  subjected  to  dynamic  pressure  loading.  A generalized, 
siagle-degree-of -freedom  model  was  utilized  because  of  its  tractability  and  because 
it  has  been  widely  accepted  as  a useful  engineering  tool.  Work  hardening  of  the 
structural  material  was  taken  into  account  to  arrive  at  a realistic  model  for  structural 
behavior.  The  analysis  considers  two  loading  conditions;  side-on  and  end-on.  Side-on 
loading  deforms  columns  and  crane  rail  beams,  while  end-on  loading  causes  rotation  and 
drop  of  the  crane  itself. 

Uncertainties  in  the  details  of  construction  of  the  building,  in  the  loads 
Induced  in  the  structure,  and  in  the  response  of  the  structure  to  these  loads  are 
estimated,  mid  propagated  throughout  the  analysis.  The  magnitude  of  the  uncertainties 
result  in  various  levels  of  confidence  in  the  success  of  an  attack.  Weapon  yield  and 
accuracy  requirements  are  also  determined. 

3.1  TARGET  DESCRIPTION 

The  target  selected  for  this  study  (Figure  3-1)  is  a typical  medium 
size  industrial  building,  as  described  in  Reference  4.  It  is  a single-story, 
single-span,  multi-bay  building  with  a heavy  steel  frame  designed  for  an  industry 
involving  heavy  equipment  and  the  manufacture  of  large  products.  The  typical  span 
length  for  this  class  of  building  is  80  ft  and  the  bay  size  is  about  40  ft.  A 50-ton, 
column-and-girder-supported  travelling  crane  is  mounted  on  beams,  bridging  the  span  of 
the  buildings.  Steel  columns,  built  up  from  rolled  steel  sections,  carry  the  roof  and 
crane  loads,  and  either  concrete  or  steel  roof  trusses  support  tha  main  roof  load 
carrier.  The  joists  on  both  the  top  and  bottom  chords  are  braced  to  prevent  lateral 
buckling. 
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Precast  reinforced  concrete  foundations  provide  seats  for  the  columns.  When 
the  column  is  Inserted  and  fastened,  this  joint  can  carry  moments  transferred  from  the 
columns.  The  foundations  extend  about  six  feet  below  the  finished  grade  to  prevent 
frost  action  during  winter.  The  crane  support  girders  are  rolled  steel  beams  or  steel 
plate  girders.  The  walls  are  supported  by  the  foundation  beams  which  bridge  the  span 
between  column  foundations.  From  the  ground  level  up  to  about  ten  feet  above  the  ground, 
tae  walls  are  made  of  brick  or  large  concrete  blocks.  Above  this  height,  the  walls  are 
made  of  light  weight  concrete  panels  or  asbestos-cement  corrugated  sheeting.  Commonly 
used  roofing  materials  are  steel  sheet,  asbestos-cement,  and  tile. 
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3.2  STRUCTURAL  RESPONSE  ANALYSIS 

When  the  building  described  in  Section  3.1  is  subjected  to  a nuclear  attack, 
several  vulnerabilities  are  possible.  However,  dynamic  pressure  and  ground  shock  are 
the  most  critical  weapon  effects  for  above-che-aurfacn  structures.  Overpressure  and 
radiation  would,  of  course,  be  lethal  to  personnel  inside  the  building,  but  are  not 
likely  to  cause  permanent  damage.  Overpressure  propagates  to  the  inside  of  a building 
of  this  type  faster  than  structural  response,  minimising  the  effect  of  differential 
overpressure  between  the  inside  and  the  outside.  At  the  overpressure  levels  of  interest, 
structural  materials  are  not  affected  by  radiation,  and  cratering,  debris,  and  thermal 
radiation  are  not  significant  threats  to  the  building  integrity.  Electromagnetic 
pulse  might  damage  electrical  equipment  and  controls,  but  they  can  probably  be  repaired 
in  a few  days,  or  means  can  be  found  to  operate  without  them. 

Dynamic  pressure  acts  on  the  side  or  end  of  the  building  like  a hurricane, 
producing  a lateral  force  proportional  to  the  exposed  area.  Ground  Bhock  produces 
oscillating  motion?.,  much  like  an  earthquake,  which  result  in  loads  dependent  on  the 
building  mass,  and  possibly  cause  differential  settlement  of  columns.  Dynamic  pressure 
is  the  most  significant  loading  mechanism  for  a building  of  this  type,  since  there  is 
a relatively  large  area  exposed  to  the  wind,  and  a relatively  small  roof  mass.  Building 
and  crane-beam  damage  resul.ing  from  dynamic  pressure  are,  therefore,  the  prime  failure 
modes  for  the  industrial  building  under  consideration  when  subjected  to  a nuclear  attack. 

If  parts  of  the  roof  or  siding  are  blown  off  by  nuclear  wind  pressure, 
emergency  repairs  could  be  made  to  permit  'resumed  operations  within  a few  days.  However, 
structural  damage  to  the  crane  or  its  supporting  structure  is  likely  to  require  materials 
which  are  in  short  supply  in  a war  situation,  and  extensive  iabrication  and  erection  of 
major  structural  members.  It  was  assumed  that  use  of  the  craue  is  essential  to  the 
industrial  process  carried  on  within  the  building,  and  that  rendering  the  crane 
Inoperative  is  a suitable  measure  of  a successful  attack. 

One  mode  of  failure  for  the  overhead  crane  may  be  lateral  displacement  of  the 
crane  rail  beams,  just  enough  to  allow  the  crane  to  fall,  while  not  necessarily  exceeding 
the  beam  yield  stress.  Another  failure  mode  may  Involve  substantial  plastic  deformation 
of  the  beams.  Based  on  these  failure  modes,  two  limiting  attack  conditions  can  be 
considered:  a low  level  attack  of  sufficient  intensity  to  produce  a lateral  beam 
deflection  Just  large  enough  to  allow  the  crane  to  drop,  provided  it  is  parked  at  the 
location  where  the  maximum  lateral  deformation  occurs;  or  a higher  level  attack  in  which 
the  maximum  plastic  deflection  of  the  beam  has  been  exceeded,  causing  the  beam  to  full, 
and  drop  the  crane  from  any  parking  position  along  the  rail  beam. 
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The  geometry  of  a typical  industrial  building  suggests  that  two  primary 
loading  orientations  should  be  considered,  ride-on  and  end-on.  For  the  side-on  caBe, 
the  detonation  point  is  located  to  the  side  of  the  building,  but  not  necessarily 
centered  between  the  columns.  This  orientation  primarily  loads  the  sides  of  the 
building,  with  a much  smaller  load  on  the  ends  of  the  building. 

In  the  side-on  loading  case,  the  siding  is  assumed  to  fail  first,  since  it  is 
not  designed  to  support  large  side  pressures  across  the  span  between  columns.  This 
leaves  the  columns  and  crane  rail  beams  exposed  to  the  dynamic  blast  pressure.  Parts 
of  the  siding,  however,  might  stay  on  the  beams  and  columns  increasing  the  exposed  area 
and,  therefore,  the  load.  Such  loadings  tend  to  move  the  rail  beams  out  from  under  the 
crane,  permitting  the  crane  to  fall. 

For  the  end-on  loading  case,  the  detonation  point  is  located  in  front  of  the 
building,  not  necessarily  centered,  so  as  to  give  the  major  loading  to  the  end  with  a 
smaller  loading  acting  on  the  sides  of  the  building.  In  this  loading  case,  blast 
pressures  act  over  the  side  of  the  crane.  Even  though  such  loadings  could  subject  the 
crane  to  pressures  which  are  essentially  uniform  along  its  length,  yaw  moments  would  be 
created  because  of  the  unaymmetrical  frontal  area  exposed  to  the  blast  wave.  Such 
moments  induced  yaw  rotation,  which  may  be  large  enough  to  rotate  the  ends  of  the  crane 
oft  the  rail  beams  allowing  the  crane  to  fall  to  the  floor. 

The  free  field  equations  used  to  determine  the  loads  are  given  in  Appendix  A. 
Uncertainties  in  the  modeling,  the  loads,  and  the  basic  variables  are  discussed  in 
Section  3.3.  The  models  synthesized  to  represent  side-on  and  end-on  loading  of  the 
building  are  discussed  in  the  following  sections,  with  details  given  in  Appendices  B to  D. 

3.2.1  Sido-on  Loading 

Side  pressure  on  the  industrial  building  will  cause  rotation  of  the  foundations, 
bending  of  the  columns  as  cantilevers,  and  lateral  bending  of  the  beams.  The  model  of 
the  building  used  in  this  study  consists  of  two  vertical  columns  attached  to  the  ground 
through  flexible  torsional  spring  support.  The  rail  beam  spanning  the  distance  between 
the  tops  of  the  columns  has  support  stiffness  between  simple  and  clamped  support 
conditions.  The  top  of  the  column  is  assumed  to  carry  a portion  of  the  roof  maBS.  When 
these  structural  elements  are  subjected  to  the  dynamic  pressure  loading,  bending  and 
rotation  are  induced  in  the  vertical  columns  and  the  horizontal  rail  beams. 

The  structural  response  is  determined  by  first  approximating  the  actual  system 


by  the  elasto-plastic,  single-degree-of-£reedom  system  model,  whoBe  governing  equations 
are  summarized  in  Appei.dix  B.  In  this  model,  the  generalized  mass,  stiffness,  and  force 
are  related  to  the  mass  and  elastic  properties  of  the  building,  and  the  dynamic  pressure 
loadings  acting  on  it,  by  the  equations  developed  in  Appendix  C. 

After  the  system  peak  response  to  a given  attack  condition  was  determined,  two 
failure  boundaries  were  computed.  The  first  fe^ure  boundary  was  based  on  the  elastic 
lateral  deformation  of  the  beam  just  sufficient  to  allow  the  crane  to  fall,  if  it  was 
located  midway  between  the  columns  where  the  maximum  lateral  displacements  occurred.  The 
second  failure  boundary  was  based  on  the  maximum  stress  in  the  beam  exceeding  the  plastic 
failure  stress,  which  would  cause  crane  damage  regardless  of  itB  location. 

Figure  3-2  shows  the  damage  distances  for  both  failure  criteria,  based  on  the 
set  of  median  analysis  parameters,  given  in  Appendix  C.  It  is  seen  that  there  could 
be  a significant  enhancement  in  lethal  radius  if  the  crane  were  parked  midway  between 
columns.  However,  in  the  vulnerability  analysis  it  was  assumed  that  the  crane  was 
parked  in  its  hardest  position,  near  a column,  and  the  plastic  failure  stress  criterion 
governed.  This  assumption  is  reasonable  for  targeting  since  the  parking  position  would 
not  be  expected  to  be  known.  < 

3.2.2  End-on  Loading 

Figure  3-3  shows  the  approximate  failure  boundary  for  the  end-on  loading  case. 
This  loading  case  makes  use  of  a simpler  model  to  determine  when  crane  drop  will  occur. 

In  the  development  of  the  model,  it  was  assumed  that  the  end  wall  will  be  destroyed, 
and  the  dynamic  pressure  will  act  on  the  face  of  the  crane  itself.  The  dynamic  pressure 
will  move  the  crane  along  the  rail  support  beams,  and  at  the  same  time  induce  a yaw 
rotation  which  will  tend  to  rotate  the  ends  of  the  crane  off  the  rail  beams  until 
it  falls  to  the  floor.  This  yaw  rotation  induced  by  the  dynamic  pressure  is  a consequence 
of  the  lack  of  symmetry  in  exposed  area  and  loading  on  the  crane.  The  failure  boundary 
in  this  case  could  be  based  on  the  actual  crane  rotation,  and  the  rotation  required  to 
cause  the  crane  to  fall  between  the  two  rail  beams  at  each  end  of  the  crane.  However, 
it  takes  only  a relatively  small  rotation  for  the  crane  to  be  free  to  fall.  Therefore, 
the  failure  boundary  was  approximated  by  assuming  the  peak  disturbance  yaw  moment  to  be 
equal  tc  the  peak  resistance  yaw  moment  produced  by  forces  at  the  crane-rail  beam  in- 
terfaces. A detailed  description  of  the  model  and  the  median  values  of  the  parameters 
used  is  given  in  Appendix  D. 
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3.3  STRUCTURAL  RESPONSE  UNCERTAINTIES 

From  the  attacker's  point  of  view,  uncertainties  leading  to  various  degrees  of 
confidence  in  a successful  attack  on  the  industrial  target  may  be  grouped  as  follows: 
uncertainties  in  weapons  effects,  uncertainties  in  intelligence  data,  and  imperfect 
modeling.  In  the  analysis  which  follows,  these  uncertainties  are  combined  with  the 
random  distribution  of  miss  distances  to  obtain  probabilities  of  successful  attacks 
at  various  confidence  levels.  To  provide  a better  comparison  with  World  War  II  damage 
from  Japan,  the  building  in  this  analysis  was  assumed  to  be  oriented  side-on  where  the 
damage  criterion  is  structural  failure  of  the  rail  beam. 

3.3.1  Uncertainty  in  Loads 

The  empirical  relationships  between  overpressure  and  range  for  a given  weapon 
size  have  been  well  established  from  weapon  teats  at  the  Nevada  Test  Site  and  the 
Pacific  Proving  Ground.  The  dynamic  pressure  relationships  are  not  so  well  defined 
empirically.  The  procedure  used  in  this  study  was  to  estimate  overpressures  from  the 
empirical  Erode  equation,  and  then  compute  the  dynamic  pressure  by  means  of  the 
Rankine-Hugonlot  equation  (See  Appendix  A).  Overpressures  were  calculated  with  a 2 -o 
uncertainty  range  of  1.2.  A 2-o  uncertainty  range  of  1.2  means  that,  if  the  overpressure 
uncertainties  are  lognormally  distributed,  about  95X  of  the  overpressure  values  are 
expected  to  fall  within  the  range  of  median/1.2  to  median  x 1.2. 

3.3.2  Uncertainties  in  Intelligence  Data 

Intelligence  data  on  targets  of  interest  is  necessarily  uncertain.  This 
category  includes  the  dimensions,  thicknesses,  and  other  properties  of  the  various 
structural  elements.  The  uncertainties  due  to  intelligence  limitations  utilized  in 
this  study  for  purposes  of  illustration  are  as  follows: 

Parameter  2o  Uncertainty  Factor 

Beam 


Moment  of  Inertia 

2.0 

Depth 

1.2 

Weight 

1.3 

Column 

Moment  of  Inertia 

1.7 

Frontal  Depth 

1.0 

Weight 

1.5 

Roof 

Mass 

2.2 
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3.3.3  Uncertainty  in  Target  Model  and  Response 

Even  if  the  building's  characteristics  were  perfectly  known,  response 
prediction*  to  a given  dynamic  pressure  would  have  large  associated  uncertainties.  The 
single-degree-of-freedom  approximation  used  here  has  been  widely  accepted  as  a useful 
engineering  model.  Predictions  based  on  this  model  provide  engineering  estimates  in  the 
linear  deformation  range.  Beyond  the  yield  point,  however,  predictions  are  less 
accurate,  since  the  behavior  of  engineering  structures  up  to  total  collapse  is  difficult 
to  predict.  Other  uncertainties  included  in  this  category  are  the  degree  of  fixity 
afforded  by  Joints,  the  drag  coefficients  associated  with  the  various  shapes,  and  the 
amount  of  wall  siding  that  might  stay  attached  to  columns  and  beams,  increasing  their 
exposed  area.  The  modeling  uncertainties  considered  in  this  study  are  sb  follows; 

Parameter  2o  Uncertainty  Factor 

Beam 


Drag  Coefficient 

1.8 

End  Fixity 

3.2 

Area  Factor 

2.2 

column 

Drag  Coefficient 

1.8 

End  Fixity 

3.2 

Area  Factor 

2.2 

Ductility 

3.0 

Work  Hardening  Parameter 

3.2 

3.3.4  Uncertainty  Propagation 

Figure  3-4  shows  the  uncertainty  in  target  hardness,  in  terms  of  the  peak 
dynamic  pressure  that  causes  failure  (Q) , which  results  from  the  uncertainties  discussed 
above.  The  data  for  this  figure  were  obtained  by  expanding  the  expression  for  dynamic 
pressure  at  failure  about  the  means  of  the  various  parameters  in  a Taylor  series,  and 
truncating  the  aeries  to  the  linear  terms.  The  resulting  expression  for  the  variance  of 
Q Involves  its  partial  derivatives  with  respect  to  each  parameter  and  the  uncertainty 
of  the  parameter.  Tne  derivatives  were  calculated  numerically,  since  there  is  no  closed 
form  expression  for  f>.  Assuming  a lognormal  distribution  for  Q,  10  and  90%  confidence 
levels  were  calculated  for  the  range  of  weapon  yields  of  interest  and  were  plotted  in 
Figure  3-4.  A 90%  attacker's  confidence  means  that  there  is  a 90%  probability  that  the 
true  (but  unknown)  hardness  is  equal  to  or  less  than  the  stated  value. 
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Also  shown  in  Figure  3-4  is  a comparison  between  these  analytically  derived 
data,  and  the  actual  damage  incurred  by  similar  buildings  due  to  the  nuclear  detonations 
In  Japan  during  World  War  II.  A recent  DNA-sponsored  study  by  Lulejian  (Reference  5) 
has  provided  a detailed  classification  and  analysis  of  the  damage  data  compiled  by 
survey  teame  at  the  detonation  sites  in  Hiroshima  and  Nagasaki. 

The  structural  damage  to  the  Japanese  buildings  was  reported  in  terms  of 
the  fraction  of  damage  to  the  walls  and  to  the  roofs  of  the  buildings  and  the  distance 
of  the  building  from  the  detonation.  For  this  comparison,  the  structural  damage 
to  the  building  walls  most  nearly  parallels  the  measure  of  building  vulnerability  used 
in  the  present  study. 

The  data  base  includes  41  single-story  heavy  steel  fra-ie  buildings,  38  in 
Nagasaki  and  only  three  in  Hiroshima.  These  buildings  are  divided  into  two  groups, 
those  with  cranes  less  than  25  tons,  and  those  with  cranes  greater  than  25  tons. 

The  single-story  buildings  with  heavy  cr-nes  comparable  to  the  industrial  target 
studied  here  existed  only  in  Nagasaki,  and  only  six  data  points  are  available. 

Assuming  a yield  of  22  KT  and  a height  of  burst  of  1640  ft  for  the  Nagasaki 
detonation,  as  recommended  by  Lulejian,  the  distance  from  the  burst  of  each  of  the 
six  applicable  data  points  was  converted  to  a dynamic  pressure  and  plotted  in  Figure  3-4. 
Each  data  point  corresponds  to  a percent  damage  to  the  walls,  as  interpreted  by  the 
survey  teams.  It  is  c icted  that  the  order  of  50%  damage  can  be  used  for  general 
comparison  with  the  analytical  results  of  Figure  3-4. 

The  10  and  90%  confidence  bounds  on  the  data  for  50%  damage  to  the  walls  of  the 
building,  calculated  using  the  maximum  likelihood  estimate  method  suggested  by  Lulejian, 
are  also  shown  in  Figure  3-4.  These  confidence  bounds  are  seen  to  be  in  general  agreement 
with  the  TRW  calculated  10  to  90%  confidence  bounds  on  damage  to  the  rc.il  beams  sufficient 
to  drop  the  crane.  However,  since  the  damage  definition  and  the  buildings  differ  from 
those  used  in  this  study,  much  significance  cannot  be  attached  to  the  comparison  without 
further  study  of  the  responses  of  the  Nagasaki  buildings. 

3.4  COMPARISON  UITH  PVN  SYSTEM 

The  physical  Vulnerability  System  (Reference  6)  gives  ratings  of  a range  of 
targets  in  terms  of  hardness  at  the  reference  yield  of  20  KT,  with  a correction  factor 
for  other  yields.  The  equivalent  Physical  Vulnerability  Numbers  (PVN)  for  the 
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industrial  target  under  consideration  were  computed  to  provide  a comparison  with  these 
data.  Curves  for  targets  with  various  sensitivities  to  yield,  so-called  K-factors,  were 
calculated  from  the  formulas  given  in  Reference  6,  and  were  matched  with  Figure  3-4. 

F*>r  the  best  estimate  curve,  the  closest  PVR  number  is  12Q9.  For  the  curves  at  10 
and  90%  confidence,  the  corresponding  PVN  numbers  are  8Q9  and  16Q9,  respectively,  where 
confidence  is  defined  as  the  estimated  likelihood  that  parameters  having  systematic 
uncertainty  combine  to  bias  the  results  to  the  extent  Indicated. 

A comparison  was  also  made  between  the  procedures  used  in  this  work  and  the 
Newmark  method  (Reference  7).  For  a yield  of  10  KT,  the  Newmark  method  gives  a dynamic 
pressure  hardness  of  2.7  psi,  which  compares  well  with  the  value  of  2.9  psl  obtained  by 
the  method  presented  in  this  report  for  the  same  parameter  values.  The  small  difference 
is  due  mainly  to  the  fact  that  work  hardening  was  included  in  the  present  work,  but  is 
neglected  in  the  Newmark  method. 

3.5  ATTACKER'S  REQUIREMENTS 

The  weapon  yields  and  accuracies  required  for  kill  probabilities  (P^)  of  0.5 

and  0.9  are  shown  in  Figures  3-5  and  3-6.  These  figures  are  based  on  the  target  hardnesses 
shown  in  Figure  3-4,  and  on  a Rayleigh-distributed  miss  distance.  The  uncertainty  bands 
shown  correspond  to  the  effect  of  the  systematic  uncertainties  on  target  hardness  and 
damage  distance. 

For  a 10-7  weapon,  for  example,  a CEP  of  1000  ft  is  needed  for  a 0.9 
probability  of  severely  damaging  the  building  (Figure  3-6).  This  is  the  best  estimate, 
or  50%  confidence  value.  To  obtain  90%  confidence  that  the  0.9  P^  will  actually  be 
achieved,  a more  accurately  guided  weapon  with  a CEP  of  750  ft  is  needed.  If  systematic 
uncertainties  could  be  reduced,  the  CEP  requirement  for  a 0.9  probability  of  kill  with 
90X  confidence  may  be  relaxed. 

A reduction  in  systematic  uncertainties  has  a similar  effect  on  yield  require- 
ments. With  a CEP  of  1000  ft,  for  example,  a 20-KT  weapon  is  required  for  a PR  of  0.9 
with  90%  confidence.  This  figure  could  be  potentially  reduced  to  nearly  10  KT  If  the 
uncertainties  could  be  substantially  reduced.  It  should  be  noted  that  if  the  yield 
is  large  enough,  CEP  requirements  are  easily  met. 
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4.0  BURIED  TARGET 

Hardened  buried  structures  design/analysis  procedures  typically  consider  only  a 
moderate  amount  of  yielding,  within  the  domain  of  small  deformations,  and  apply  simplify- 
ing approximations  which  nearly  always  add  conservatism.  This  background  has  generally 
led  to  significant  underestimation  of  the  overpressure  (for  a given  yield)  required  to 
collapse  such  a structure.  In  fact,  no  fully  buried  concrete  structure  has  collapsed 
when  exposed  to  nuclear  airblast  loading;  with  a number  of  structures  retaining  their 
protective  capability  at  overpressure  higher  than  their  predicted  collapse  load. 

For  targeting  considerations  the  loading  required  to  damage  the  structure  beyond 
the  point  of  its  subsequent  usefulness  is  of  primary  interest.  In  such  cases,  proper 
modeling  of  the  large  deformation  response  is  of  paramount  importance  and  use  of  conserva- 
tisms which  have  been  applied  to  the  design  process  needs  to  be  reconsidered. 

In  analyzing  the  response  of  a rectangular  bunker  the  typical  approach  has  been 
to  model  the  roof  as  a single-degree-of-freedom  slab  with  pinned  or  fixed  supports.  When 
the  roof  deflection  reaches  a specified  factor  (ductility  ratio)  times  an  elastic 
deflection,  the  structure  is  considered  to  have  failed.  However,  the  inability  of  the 
technical  community  to  predict  structural  failure  loads  has  led  to  the  identification  of 
several  key  issues  regarding  this  approach,  including  (1)  Is  a single-degree-of-freedom 
model  adequate?,  (2)  If  so,  what  end  conditions  should  be  used?,  (3)  How  is  the  response 
affected  by  loading  on  the  bunker  sides?,  (4)  How  does  structural  damping  influence 
the  response?,  (5)  What  is  the  effect  of  rigid  body  motions  of  the  whole  bunker?, 

(6)  What  is  the  large  deflection  force-deformation  behavior  of  the  structure?,  and 
finally,  (7)  What  is  the  failure  criterion? 

Since  the  analysis  uncertainties  have  a major  impact  on  the  targeting  un- 
certainties, some  effort  was  placed  in  evaluating  the  analysis  model  to  use  in  the 
study.  It  was  determined,  based  on  comparison  with  more  detailed  finite  element 
calculational  results,  that  a single-degree-of-freedom  model  would  be  adequate.  However, 
tho  model  used  in  this  study  considered  the  primary  response  mode  of  the  complete 
structure  and  first-order  large  deflection  force-deformation  structural  behavior,  which 
includes  effects  of  in-plane  stresses.  Failure  criterion  is  based  on  the  structure 
achieving  a specified  ductility  or  deflection.  Rigid  body  motions  of  the  bunker  have 
not  been  included  directly  in  this  study.  However,  a straightforward  modification  of 
this  analysis  is  all  that  is  required  to  incorporate  inertial  loading  effects  caused 
by  these  motions. 
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In  order  to  evaluate  the  elastic  response  portion  of  this  analysis  model  the 
bunker  configuration  studied  is  identical  to  that  used  in  the  Reference  8 finite  element 
analysis  performed  by  the  Waterways  Experiment  Station  (WES).  Since  little  data  are 
available  to  properly  evaluate  the  large  deformation  portion  of  the  response,  this  study 
primarily  addresses  the  sensitivity  of  the  response  to  possible  variation  in  the  primary 
model  parameters.  Estimates  of  the  uncertainties  in  the  model  parameters  and  in  the  model 
itself,  alorg  with  the  sensitivity  coefficients, determine  the  target  hardness  uncertainty. 

4.1  TARGET  INSCRIPTION 

A box-type  reinforced  concrete  bunker  structure  is  assumed  to  be  sufficiently 
long  co  that  the  response  can  be  considered  two  dimensional.  Obviously,  short  structures 
with  the  same  cross  section  will  be  stronger.  Figure  4-1  shows  the  40-ft  wide  by  25-ft 
high  bunker  to  be  fully  buried  under  5-ft  of  overburden.  The  5-ft  thick  roof  slab  has 
2.0%  and  1.5%  reinforcement  at  the  slab  bottom  and  top,  respectively;  while  the  4-ft  thick 
walls  have  1.0%  compression  and  tension  steel.  Other  properties  of  the  bunker  structure 
and  surrounding  backfill  are  shown  in  the  figure.  The  site  medium  ia  dry  soil  over  shale 
with  the  depth  to  the  soft  rock  located  approximately  at  the  bunker  foundation  level. 

4.2  STRUCTURAL  RESPONSE  ANALYSIS 

The  bunker  is  considered  to  be  subjected  to  effects  caused  by  the  overpressure 
loading  of  a nuclear  weapon  detonation.  The  side  loading  on  the  walls  is  based  on  the 
stress  attenuation  associated  with  the  overhead  airblast  loading.  Thus,  in  this  analysis, 
direct-induced  and  upstream-airblast-induced  effects  have  been  neglected.  In  layered 
stiffer  geologies  these  effects  could  be  quite  significant  and  the  possibility  of  the  wall 
being  more  vulnerable  than  the  roof  should  be  considered. 

If  the  primary  response  mode  of  the  bunker  is  known,  then  "simplified"  analysis 
procedures  can  be  applied.  These  procedures,  which  assume  the  bunker  to  respond  only  in 
this  primary  mode,  would  be  limited  by  (1)  how  closely  the  actual  response  corresponds  to 
the  assumed  mode  and  (2)  knowledge  of  the  complete  load-deformation  relationships.  The 
roof  of  the  bunker  considered  in  this  study  is  expected  to  be  the  most  critical  vulner- 
ability, which  is  consistent  with  the  assumed  response  mode.  For  conditions  where  the 
wall  response  is  most  critical  a different  mode  would  be  considered. 

During  the  initial  phase  of  the  bunker  response  the  structure  behaves  elasti- 
cally, followed  by  various  degrees  of  plastic  deformation  and  finally  collapse.  The 
following  subsections  address  each  of  the  response  phases. 
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The  bunker  is  modeled  as  four  rigidly  connected  one-way  slabs  whose  motion  are 
primarily  due  to  bending  effects,  and  therefore,  only  the  midplnne  motion  of  these  slabs 
is  required  to  define  the  complete  bunker  response.  The  response  of  the  ith  slab  is  then 
approximated  by  (Figure  4-2) 

yi  " ^i^P  (4-1) 

where  ^(positive  Inwards)  is  the  normalized  mode  shape,  n is  the  single-degree-of-freedom 
response  amplitude  and  Is  a coordinate  along  the  midplane  of  the  ith  slab  (roof,  i ■ 1; 
walls,  i ■ 2 and  3;  floor,  i ■ 4).  This  results  in  the  undamped  equations  of  motion 
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and  p^,  D^,  and  o^  are  the  mass  density  per  unit  midplane  area,  flexurul  rigidity, 
thickness  and  transverse  applied  stress,  respectively,  of  the  ith  slab  und  E and  v ure 
Young's  modulus  and  Poisson's  ratio,  respectively,  of  the  concrete.  Substituting  the 
Equation  (4- Dapproximation  into  the  equations  of  motion,  multiplying  by  ^ , integrating 
with  respect  to  x^  and  summing  over  the  four  slabs  yields 
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(4-4) 

(4-5) 

(4-6) 


are  the  generalized  mass  and  frequency  relationships. 

Since  this  analysis  considers  a bunker  configuration  for  which  the  roof  vulner- 
ability is  most  critical,  the  elastic  mode  shape  (Table  4-1)  was  based  on  the  deflected 
shape  for  a static  uniform  pressure  (w)  loading  on  the  bunker  roof.  The  moments  and 
at  the  intersections  between  the  walls  and  roof  and  walls  and  floor,  respectively,  are 
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4—2.  Structural  Response  Analysis 


Table  4-1  Bunker  Modal  Equations 
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determined  from  condition*  of  slope  continuity  at  thu  corners,  and  the  equivalent  slab 
length*,  hj,  are  meusurud  from  the  midpoint*  of  thu  corner*. 

In  proceeding  with  thu  unuly*i*  thu  following  simplifying  approximation*  are 
made  regarding  the  loading!  (1)  thu  ma**  of  aoil  overburden  ia  considered  part  of  thj 


roof  which  is  loaded  uniformly  by  u time  varying  overpressure,  (2)  the  sidu  louding  Is  also 
uniform  and  time  varying  and  is  symmetrically  applied  on  both  sidu*  of  thu  bunkur  and  (3) 
the  effect  of  thu  bottom  louding  1*  neglected.  Then  liquation  (4-4)  becomes 


H + 2tiuf)+  | ^ 


(4-7) 


where  modal  dumping  ha*  burnt  added  with  t!  being  the  percent  critical  damping  and  thu 
gcnuralUod  mass  and  frequency  relationships  are  given  in  Tuble  4-1  (iS  being  a tracer 
purumutur  which  is  equal  to  1.0  for  thu  elastic  portions  of  the  response  and  y and  y 
are  the  weight  densities  for  concrete  and  soil  overburden,  respectively). 

liquation  (4-7)  is  directly  integrated  using  the  Urodo  Ill'll  overpressure  function 
(Kefurencc  0)  and  un  expression  for  laterul  loud  based  on  thu  stresu  attenuation  in  a 
bilinear  medium  (linear  louding  waveapeud,  e(),  and  unloading  waveapecd,  Cj)  for  a single 
exponential  approslmat ion  (p  exp  - (t/tuo))to  the  surface  pressure  (Kci erence  lli , via., 
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and  p Is  the  peak  surface  overpressure,  is  the  initial  slope  time  intercept,  k Is  the 
laterul  stress  coefficient  and  Ji  t*  the  mid-duptli  of  the  structure, 

A purely  ulastie  analysis  of  the  baseline  hunker  configuration  was  performed 
to  compare  results  of  this  simplified  analysis  with  the  more  detailed  WPS  ealeulat ional 
results.  This  comparison,  which  is  shown  in  figure  4-3,  indicates  the  peak  deflection  tn 
the  WK8  results  to  be  slightly  larger  than  in  the  simplified  analysis  (.‘.b  tn.  as  compared 
to  2.7  in.).  The  waveform  comparison  is  good  for  the  I list  response  cycle;  with  the  Wliii 
calculation  exhibiting  larger  equivalent  damping  at  later  cycles,  a*  might  be  expected,  in 

both  cases  the  response  oscillate*  about  the  qua*l-*tatic  solution  (static  detraction 
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caused  by  the  overpressure  at  a given  time) 


A comparison  between  the  simplified  analysis  and  WES  results  for  their  stiffer 
backfill  case  does  not  show  good  agreement;  probably  caused  by  the  approximations  used  in 
calculating  the  side  loading  (e.g.  the  side  loading  does  not  act  until  the  backfill  stress 
wave  gets  to  the  midpoint  of  the  wall).  Thus  modeling  of  this  loading  needs  further  study. 
However,  the  baseline  configuration  of  this  study  considers  the  soft  backfill  whose  load- 
ing does  not  have  a significant  effect  on  the  bunker  roof  response. 

4.2.2  Small  Deformation  Plastic  Response  Phase 

The  elastic  analysis  is  considered  to  be  applicable  until  the  fully  plastic 
moment  develops  at  the  top  or  bottom  corners  or  at  the  midpoint  of  the  roof.  This  dis- 
cussion assumes  that  the  plastic  hinge  initially  develops  at  the  top  corners  (the  case  for 
the  baseline  buuker  configuration),  however,  similar  analysis  methodology  would  be  applic- 
able if  yielding  initiated  at  the  center. 

Yielding  is  governed  by  the  moment-thrust  diagram  (Reference  11)  shown  in 
Figure  (4-4).  The  diagram  is  a linear  fit  between  (1)  the  compression  capability  of  the 
slab  with  zero  moment,  Pq,  (2)  the  balance  point  M^,  Pfa,  corresponding  to  the  maximum 
moment,  (3)  the  moment  capability  in  pure  bending,  Mp  and  (4)  the  tension  capability  of 
the  slab.  For  compactness  the  equations  governing  the  points  on  the  yield  surface  are 
given  in  the  figure. 

It  is  not  known  apriori  whether  yielding  at  the  roof  corners  will  occur  within 
the  roof  or  wall  slabs  and  therefore  both  possibilities  must  be  checked.  The  moment  is 
determined  from  the  curvature  at  the  corner  as 


Ma  - <l-k1)n(t)  (4-9) 

L1 

while  the  axial  load  in  the  roof,  determined  from  the  shear  (dM/dx)  in  the  wall  just  below 
the  corner,  is  given  by 


A1 


1 - 


kl  + 


6k2k3 


n(t) 


l_k d 

2 2 


(4-10) 


and  similarly  for  the  axial  loads  in  the  walls 


(4-11) 
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The  ultimate  moment  is  reached  when  M^,  or  MA,  P^reach  the  yield  surface. 

When  the  yield  hinge  develops  the  stiffness  and  primary  frequency  of  the 

structure  decrease  and  the  mode  shape  must  be  modified.  The  load,  w,  used  as  a basis  for 

calculating  the  revised  mode  shape  is  taken  as  the  average  of  the  static  loads  required 

Y 

to  (1)  develop  the  yield  moment,  M.,  at  the  corners  and  (2)  develop  the  yield  moment, 
y Y A 

Mc,  at  the  c‘>r't*-r  with  applied  at  the  coiners.  The  mode  shape  is  then  proportional 
to  the  static  deflection  caused  by  this  load  and  by  the  corner  moments  M^.  Thus,  the 
equation  £01  in  Table  4-1  is  modified  to  be 


(4-12) 


for  purposes  of  calculating  the  mode  shape  <j)^,  and  the  corresponding  equation  of  motion 
is 


H + 2f5wh 


. 2 

+ « n 


(4-13) 


where  the  generalized  mass  (M^)  and  frequency  (u)  are  as  defined  in  Table  4-1,  with  the 
tracer  parameter  6 equal  to  zero.  Continuity  of  displacement  and  velocity  at  the  center 
of  the  roof  is  the  basis  for  determining  initial  conditions  for  the  plastic  region.  This 
small  deformation  plastic  response  formulation  is  considered  to  be  applicable  until  the 
fully  plastic  moment  is  developed  at  the  center  of  the  roof. 


Subsequent  response  can  result  in  large  deformation  since  three  yield  hinges 
have  developed.  However,  resistance  to  this  motion  is  still  provided  by  the  yield  moments 
and  by  the  in-plane  tensile  forces  which  develcp  as  the  roof  deformation  becomes  large. 

4.2.3  Large  Deformation  Plastic  Response  Phase 

Once  three  yield  hinges  develop  in  the  roof  then  bending  effects  between  hinges 
will  be  neglected.  The  equation  of  motion  of  the  link  between  two  hinges  is  given  by 


with 


0 


(4-14) 
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(4-15) 
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T " ES  (Pi  + PpHj^U  + 0^)  - lj 


(4-16) 


is  the  tension  in  the  roof  slab,  where  Ec  is  the  steel  modulus,  and  0 is  the  rotation  for 

a Y Y 

the  link.  During  this  phase  of  the  response  the  yield  moments  and  are  reduced 
according  to  the  moment  - thrust  criterion  of  Figure  (4-4).  When  the  maximum  tension, 

Tq,  is  reached  the  moments  go  to  zero  and  the  tensile  load  thereafter  is  given  by 


To  " fy(Pl  + Pi>Hi 


(4-17) 


where  f is  the  yield  stress  of  the  steel. 

y 

4.2.4  Structural  Capability 

Figure  (4-5)  shows  the  response  of  the  baseline  bunker  configuration  for  a 
l-MT/300-psi  overpressure  loading.  Initial  yielding  occurs  at  the  top  of  the  wall  at 
approximately  0.013  secs.  If  the  subsequent  response  followed  the  yield  surface  of 
Figure  (4-4)  the  wall  moment  would  continue  to  increase  since  the  axial  load  increases, 
from  a state  below  the  balance  point.  However,  whether  this  additional  capability  is 
achievable  after  yielding  occurs  is  questionable  and  the  yield  moment  has  been  taken  to 
be  constant  when  the  post-yielding  axial  compression  increases  (below  the  balance  point). 
The  yield  surface  is  followed  when  the  axial  load  decreases,  as  shown  in  Figure  (4-4), 

The  moment  in  the  center  of  the  roof  continues  to  increase  in  the  plastic 
region  until  the  yield  criterion  for  the  roof  is  satisfied.  (The  initial  small  moment 
decrease  in  this  region  is  an  inaccuracy  in  the  model  caused  by  the  change  in  mode  shape 
at  the  beginning  of  the  first  plastic  region.)  At  0.016  secs  three  yield  hinges  have 
developed  in  the  roof  and  large  deformations  can  occur  which  cause  tension  in  the  roof 
slab.  The  residual  roof  moment  capability  associated  with  this  tension  is  indicated  in 
Figure  (4-5).  A maximum  roof  deflection  of  8.3  inches  occurs  at  time  of  0.084  secs, 
which  corresponds  to  a ductility  of  5.5,  with  respect  to  the  three  yield  hinge  deflection 
and  a ductility  of  8.3  with  respect  to  the  elastic  deflection. 

The  bunker  response  to  30  hi'/ 300  psi  and  10  KT/3000  psi  overpressure  loadings 
i3  shown  in  Figures  (4-6)  and  (4-7),  respectively.  In  the  latter  case,  the  large 
roof  tension  and  correspondingly  large  shear  stresses  at  the  top  of  the  wall  would 
actually  tend  to  be  relieved  by  cantilever  yielding  at  the  bottom  of  the  wall 
(currently  not  included  in  the  model).  However,  neglecting  this  effect  would 
not  be  expected  to  have  a significant  effect  on  the  results  since  (1)  the  moment 
caused  by  the  roof  tension  was  relatively  small  and  (2)  the  shear  loading  was  not 
considered  to  reduce  the  moment  capability  at  the  top  of  the  wall  (showo  to  be 
constant  under  large  deformation) . 
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Figure  4-6.  Bunker  Response  for  10  KX,  300  psi  Overpressure 


£ Center 


Bunker  Response  for  10  KX,  3000  psi  Overpressure 


A more  general  bunker  model  should  include  combined  large  deformation  behavior  of  the 
roof  and  walls  and  yield  conditions  which  include  the  effects  of  shear  loading  in 
addition  to  moment  and  thrust. 

Figure  (4-8)  shows  the  damage  function  corresponding  to  a baseline  ductility 
of  10  (with  respect  to  three  yield  hinge  deformation)  and  to  a center  deflection  of  10% 
of  the  roof  span  (3.6  ft  or  p ■ 28.5).  If  the  larger  deformations  were  achievable,  a 
significant  increase  in  hardness  over  the  y of  10  case  would  result.,  especially  at  the 
lower  yields  (factor  of  2.7  for  10  KT  compared  with  1.6  for  1 MT).  For  purposes  of  further 
comparison  damage  functions  are  shown  for  a roof  slab  with  simply  supported  and  with  fixed 
end  conditions;  both  for  ductilities  of  10.  The  "lower  bound"  simply  supported  condition 
reflects  hardnesses  lower,  by  a factor  of  1.7  at  10  KT  and  a factor  of  1.5  at  1 MT,  than 
determined  from  the  analysis  model  used,  with  the  corresponding  factors  for  the  "upper 
bound"  fixed  end  conditions  being  1.5  at  10  KT  and  at  1 MT.  To  investigate  effects  of 
wall  yielding  reducing  the  roof  tension,  additional  analysis  considered  the  maximum  tension 
to  be  zero  with  the  minimum  moment  equal  to  the  roof  capability  under  pure  bending  (M^). 
These  assumptions  resulted  in  a factor  of  1,2  increase  in  hardness  for  1 MT  and  1.4 
increase  for  10  KT. 

4.3  STRUCTURAL  RESPONSE  UNCERTAINTIES 

Uncertainties  in  the  bunker  response  are  associated  with  uncertainties  in  (1) 
the  structure-media  interaction  (SMI)  response  analysis,  (2)  the  structure  configuration, 
(3)  the  material  properties  and  (4)  the  free  field  weapon  effect  environments.  Since  a 
generic  bunker  design  has  been  considered,  uncertainty  sources  associated  with  intelligence 
data  have  not  been  considered.  However,  even  for  a high  degree  of  intelligence  data 
(e.g.,  structural  drawings,  material  specifications  and  average  engineering  site  pro- 
perties are  known),  uncertainties  still  exist  in  such  parameters  as  structural  material 
properties,  since  contractors  typically  use  materials  which  have  properties  better  than 
specification.  For  example,  statistical  analysis  of  two  sets  of  test  data  for  40,000  psi 
steel  reinforcing  bars,  summarized  in  Reference  (12),  shows  the  actual  strength  to  be  1.19 
and  1.27  times  the  nominal,  respectively.  Targeting  vulnerability  analysis  can  account 
for  material  properties  which  are  better  than  design  specification  requirements. 

Determination  of  the  bunker  hardness  uncertainties  is  based  on  the  methodology 
in  Section  2.2.  Sensitivities  of  the  hardness  to  variation  in  the  analysis  parameters  are 
determined  as  a function  of  yield  from  a series  of  response  calculations,  for  specific 
yields,  each  of  which  consider  variation  of  an  individual  parameter.  Table  4-2  summarizes 
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Tabl*  4-2  Systematic  Bunker  Uncertainties 


Category 

Parameter 

Structure 

Media 

Interaction 
Kesponse 
vtut  lysis 

Mode  Shape, 

Etptlv,  Hoc' T Span,  1.^ 
Failure  Criteria,  \- 

Heading  Yield,  M 

Tens  Hi  Yield,  T 

(2\  ° 
Damping,  |-  *■' 

Ana lv si  a Model 

Subtotal 

Structure 

Configuration 

Steel  Area,  p 4 p* 
Fffectivo  Hoof  Depth, 
Overburden  Depth,  11,, 

Subtotal 

Materia  1 

1'ropcr t ion 

Steel  Yield,  f 

V 

Concrete  Module"-.  K 

c 

backfill  Strewn 

Speed,  C,(1> 

Subtotal 

Weapon 
l’f  f out a 

Teak  OvrrpreHHure,  p 

Initial  Slope 

Intercept,  t 

Subtotal 

Sensitivity 

-jLiiiutl 

<Kln  x„) 


10  KT 


2a 

*i 

Uncert. 

Factor 


2 r Target 
Hurdnesfi 
Uncertainty 


10  KT 


Tot  n 1 


the  key  parameter  sensitivities  l'o*-  10  KT  and  1 MX  yield**,  the  parameter  2<r  uncertainty 
factors  and  target  Uurdnuuu  uncertainty  component,  corresponding  to  both  yields,  associated 
with  euch  parameter  and  with  the  analysis  model.  The  basis  for  estimating  the  2 o 
uncertainty  factors  in  the  unuiysis  parameters  is  summurixod  in  Tublu  4-3. 

Given  reasonably  good  intelligence  information  on  the  structure,  it  is  suen  tliut 
the  analysis  parameter  uncertainties,  particularly  thu  failure  criteria,  huvu  the  most 
significant  effect  on  turget  hardness  uncertainty.  Thu  unuiysis  model  uncertainty 
considers  that  even  if  ail  model  parameters  are  known  udequutuly,  the  unuiysis  may  not 
predict  thu  responsu  correctly.  This  uncertainty  would  be  strongly  dependent  on  thu 
bunker  being  analysed,  i.u.,  if  it  were  resting  on  soft  soil  such  that  rigid  body 
motions  wure  importa.rt  or  if  thu  side  iouding  had  a significant  effect  on  thu  response, 
thu  unuiysis  model  uncertainties  would  bu  iurgur  than  for  thu  bunker  considered  in  this 
study.  Model  uncertainties  huruin  are  considered  to  be  largest  where  modeling  of  thu 
structure  yielding  is  most  critical  (Low  yields  wuupons)  but  uru  estimated  to  bu  nut 
us  lurgu  us  the  uffucts  of  failure  critoriu  uncertainties.  Additional  teat  dutu  for  the 
lurgu  deformation  response  of  model  bunkers  is  required  to  better  dufino  thu  failure 
crituriu  und  vuiidutu  unuiysis,  especially  tor  the  low  yiuid  woupon  loading.  Current 
uncertainties  may  have  an  important  impact  on  weapon  allocation. 

Dumping  wus  not  considered  in  the  baseline  analysis  (although  it  may  have  been 
appropriate  to  consider  the  order  of  20%  damping)  und  it  is  seen  that  the  affect  of  dumping 
or  its  uncertainties  is  not  vary  significant.  The  target  hurdnsss  uncertainties  shown  in 
Table  4-2  ure  based  on  un  analyeie  which  considered  20%  dumping  within  the  smull  deforma- 
tion portion  of  the  rusponse.  Rather  large  Hardness  increases  are  obluiuod  if  thu  sumu 
dumping  coefficient  is  usud  in  both  the  elastic  and  plusllc  region.  However,  within 
thu  plustic  region,  thu  dumping  forces  should  net  Imply  u total  resistance  which  Is 
signif icuntiy  beyond  the  structure  capability. 

WES  analysis  results  indicate  that  utiffur  backfill  modiu  will  reduce  the 
roof  response.  It  Is  estimated  that  this  effect  would  not  bu  greater  than  huving  kite 
roof  act  us  a fixed-end  beam.  Tlius  thu  ul'fect  of  backfill  uncertainty  on  turget  hurdnuss 
is  considered  to  be  the  same  as  the  effect  of  mode  shape  uncertainty,  in  calculating 
totui  target  hurdness  uncertainty  this  effect  is  counted  only  once,  which  iH  equi- 
valent. to  huving  a correlation  coefficient  of  -0,5  between  the  target  hardness 
uncertainties  corros ponding  to  these  parameters . Although  usud  in  this  study,  a 
lognormal  distribution  is  not  un  accurute  representation  of  the  backfill  roluted 
uncertainty  since  softening  the  backfill  does  not  influence  the  response  nearly  as 
much  us  stiffening  does.  However,  the  approximation  will  not  affect  the  totui 
uncertainty. 


45 


n iti  ,.ri^KW 


Tibi*  4-3  Uiaid  For  P»r*m*t«r  Uncertainty  Eatimatni 


1 


Catagory 


St  rueturu 
Mod  la 

Interact  Ion 

KciipOUNU 

Analyala 


Parana ter 

x. 


Structure 
Coni' inanition 


Material 
Proport  tea 


Weapon 
If  feet  a 


Mode  Shape,  if 


lojulv,  Hoof  Span,  1 


'1 


Failure  Criteria,  p 


Hand  1 tut  Vi  eld,  M. 


‘Iona  He  Vie  Id,  T 


Damp  Inn,  II 


Aaalyalu  Model 


Steel  Area,  p A p* 


Kfl'eetlvo  Hoof  Depth,  II 


Overburden  Depth,  1L 

o 


Steel  Yield,  f 

V 


Concrete  ModuluH,  K 


lluekt'lll  StroHH  Speed,  C, 


Peak  Overproduce,  p 


Initial  Slope  Intercept , t 


Uncertainty  Factor  bauia 


Fixed-end  and  aimply-Bupported 
reel'  beaut  anal  y a la  re  unit  a 


Width  it l'  walla 


Difference  betwuun  u M 10  and 
0,10  Lj  roof  deflection  reanlta 


Heforenao  (8) 


Concrete  tension  capability 


Kffeot  of  2 OX  damping 


lintimatod  baaed  on  failure 
criteria  related  unccrvutntiuH 


linn  inert  iny,  eatimate 


Variation  In  at  eel  oover 


Dew inn  depth  of  40  lu.|  nominal 
of  f>0  tn,t  maximum  of  *)0  In. 
fitr  the  median  value 


Reference  (S) 


Reference  (H)i  K(i  •v  \f 1^', 


Would  not  have  greater  effect 
than  fixed-end  slab 


Reference  (')) 


Same  an  p 


Tho  total  target  hardneuu  uncertainty  lu  eonulderubly  larger  for  10  KT  wuuponu 
(2#  factor  of  4.0)  than  for  1 MT  wuuponu  (2*  factor  of  2.3).  Thin  iu  primarily  uuuociutud 
with  tho  fuel  that  at  high  yield  tho  ovorpruuuuro  required  to  fail  tho  utructuru  iu 
not  uuuriy  uu  dupondont  on  the  ductiiity  uu  iu  tho  cuuo  for  low  yiold.  On  tho  buulu 
of  thouo  rcuuitu  damage  function  curvuu  corrouponding  to  10%,  50%  (bout  outimutc) 
and  00%  confidoucu  are  uhown  In  figure  4-9,  Tho  30%  confidence  curve  iu  calculated 
from  tho  modiun  vuiuou  of  ail  unuiyuiu  parumotoru  and,  baued  on  related  anulyuiH,  iu 
only  very  ulightly  different  from  tho  uotunl  median  hardneuu. 

4.4  COMPARISON  WITH  l'VN  SYSTEM 

for  tho  buut  outlmutu  curve,  the  cioeout  l'VN  (Reference  b)  lu  3/1’’/.  for  tho 
10  and  '»0%  confidence  curvou,  tho  corrouponding  numberu  are  and  411*7,  rouped  fvely, 

4 . 5 ATTACKER ' S REQU iREMENTS 

Attacker  roquiromontu  generally  conuider  prohub ilitiuu  of  kill,  l*R,  of  much 
greatur  than  50%.  Rundom  vuriutionu  will  have  uomu  influence  on  thouo  roquiromontu 
uinco  tho  targeting  may  need  to  conuider  neur-worut-uite  conditionu.  However,  thuuo 
uto  generally  negligible  compared  to  the  random  vuriutionu  introduced  by  the  CEP 
diut  ributiou. 

Table  4-4  givuu  eutimutou  of  random  variation  effoctu  for  both  thu  10  KT  and 
1 MT  yloldu,  bailed  on  thouo  dutu  and  on  tho  ayatomutic  uncertainty,  the  weapon  yiuld 
and  CEP  accuracy  required  for  P^,  of  0.5  and  0,9  wuu  calculated.  Tho  t'ouultu  are 
uhown  in  figurou  4-10  and  4-11  for  confidence  levolu  of  10,  30  und  90%.  Table  4-4 
glveu  eutimutou  of  random  variation  olfoctu  for  both  the  10  KT  und  l MT  yieldu. 

However,  uh  diucuuuod  in  Section  2,  thoue  effoctu  do  not  tiiguif icuntly  influence 


Category 


Structure 
Conf Igurntlnn 


Materia  1 
Propert  lea 


Weapon 
13  f f e e t k 


Parameter 

x 

Sensitivity 
d(ln  p) 

2o  1 

Xi 

2o  Target 
Hardness 
Uncertainty 

d(Tn 

“V 

Uncert. 

i 

10  KT 

1 MT 

Factor 

10  KT 

1 MT 

Steel  Area,  p and  p' 

1.45 

1 .00 

1.0 

1.0 

1.0 

Effective  Roof  Depth, 

2.0 

1 .8 

1.04 

i . 08 

1.07 

Overburden  Depth,  11^ 

0.22 

0.11 

1.5 

3.09 

1.04 

Subtotal 

1.3 

3.1 

Steel  Yield,  f 

y 

1.34 

1.00 

3.2 

1.3 

1 .2 

Concrete  Modules,  K 

0,40 

0.20 

1.2 

1 . 08 

1.04 

Haekfill  Stress  Speed, 

C 

0 

1.6 

1.5 

Subtotal 

3.7 

1.6 

Peak  Overpressure,  p, 

0.73 

0.97 

1.1 

1.07 

1.10 

Initial  Slope  Intercept, 

t 

r»' 

1.52 

0.77 

i 

3 . 1 

3.16 

3 .08 

Subtotal 

■a 

1.1 

5.0  CONCLUSIONS 


This  study  demonstrates  the  use  of  an  analytical  methodology  for  the  evaluation 
of  the  uncertainties  associated  with  the  prediction  of  target  vulnerability.  Two  types 
of  uncertainties  are  properly  distinguished,  namely:  random  variations,  which  represent 
actual  nhysical  variations;  and  uncertainties  which  represent  potential  bias  between 
the  real  world  and  the  perceived  world  of  the  analytical  models  and  model  parameters. 

Random  variations  affect  the  damage  probability  of  a given  target:  (or  the 
expected  number  of  damaged  targets  for  attacks  on  similar  targets)  and  are  treated 
adequately  by  the  Physical  Vulnerability  System.  However,  the  more  significant  systematic 
uncertainties  are  not  treated  adequately  by  this  system.  Where  high  confidence  of  high 
target  damage  probability  is  required,  the  targeting  requirements  for  most  cases  are 
governed  by  the  systematic  uncertainties. 

The  targeting  uncertainties  for  the  industrial  building  studied  are  largely 
associated  with  the  unique  construction  of  a heavy  crane  building.  These  uncertainties 
could  be  reduced  by  experiments  designed  to  provide  force-versus-def lection  data  for 
structures  similar  to  the  50-T  crane  building  analyzed.  However,  the  cost-to-benef it 
ratio  for  such  experiments  would  be  prohibitive  unless  correlations  could  be  established 
between  the  representative  structures  tested  and  a number  of  industrial  targets. 

The  targeting  uncertainties  for  the  buried  bunker  analyzed  are  largely 
associated  with  the  fact  that  the  analytical  procedures  have  not  been  validated.  The 
results  of  this  study  indicate  that  simplified  analytical  models  which  include  the 
primary  response  mechanisms,  such  as  the  effect  of  wall  flexibility  and  post-yielding 
resistance,  are  adequate.  An  advantage  of  a mechanistic  model  over  a finite  element 
analysis  is  that  it  is  not  restricted  to  a continuum  description  of  the  structure 
analyzed. 

Additional  test  data  are  needed  in  several  areas  to  reduce  the  uncertainties 
in  the  present  analysis.  Dynamic  testing  should  consider  representative  construction 
techniques,  backfill  placement,  and  site  media  where  possible.  Initial  testing  should 
consider  the  smallest  size  models  for  which  these  conditions  could  be  satisfied.  The 
airblast  loading  should  be  simulated  for  times  corresponding  to  structural  failure. 

Both  high  and  low  yield  testing  is  desirable,  but  the  high  pressure,  low  yield  end  of 
the  hardness  curve  should  be  emphasized  since  the  uncertainties  are  largest.  Testing 
at  the  high  yield  end  provides  baseline  hardness  data. 
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Both  the  simplified  analysis  model  and  appropriate  finite  element  analysis 
should  be  vised  for  pre-test  production  and  post-test  evaluation.  To  support  this 
effort,  several  features  could  be  added  to  the  simplified  model,  including  (1)  Improved 
treatment  of  the  bunker  side  loading,  (2)  rigid  body  motion  of  the  bunker,  and  (3) 
effects  of  soil  arching.  The  specific  modeling  of  the  soil-media-interaction  phe- 
nomenology needs  to  be  verified  by  comparison  with  results  of  finite  element  analysis 
and  test  data. 


6.0  RECOMMENDATIONS 


Based  upon  the  results  of  this  study,  the  following  general  recommendations  are 

made: 

1)  Error  analyses  of  target  hardness  estimates  should  be  associated 
with  each  analytically  derivud  target  hardness  estimate  so  that 
the  potential  users  of  such  information  cun  evaluate  the 
associated  risks. 

2)  The  Oy  term  of  the  Physicul  Vulnerability  System  should  not  be 

used  to  include  systematic  errors  of  anulyses.  The  o term  very 
often  has  little  effect  upon  the  probability  of  damage.  Hence 
including  analytical  uncertainty  understates  the  risk  of  significant 
errors  in  estimated  dumugc  probabilities. 

3)  Additional  testing  of  scale  model  bunkers  should  establish  un 
empirical  basis  for  predicting  large  deformation  and  establish 
failure  criteria  of  this  class  of  structure. 

4)  Simple  analytical  formulations  for  buried  bunker  response  should 
be  validated  with  test  data. 

5)  Sensitivity  analyses  should  be  conducted  to  provide  a dntu  buse  to 
support  weapon  effectiveness  evaluation  and  research  planning. 
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APPENDIX  A FREE  FIELD  EQUATIONS 


A.l  Peak  Overpressure 

Equation  (3-3)  of  Reference  11  expresses  the  peak  overpressure  p 
in  psi  in  terms  of  the  yield  W expressed  in  MT  and  radial  distance 
r expressed  in  kft  created  by  a surface  burst  by  the  equation 


P 


so 


3300W  . 192  fT 

r3  r Vr 


(Al) 


A. 2 Peak  Dynamic  Pressure 


Equation  (3-10)  of  Reference  11  expresses  the  peak  dynamic  pressure 
q (psi)  in  terms  of  the  peak  overpressure  p (psi)  and  ambient 

nO  r 3 O 

pressure  pQ  (psi)  by  the  equation 


q 


o 


2.5  p 


so 


P + 7 P 
*so  *o 


(A2) 


A.  J Time  Variation  of  Dynamic  Pressure 

The  peak  dynamic  pressure  occurs  shortly  after  the  arrival  of  the 
blast  pressure  and  subsequently  decays  approximately  exponentially  with 
time.  It  is  assumed  here  the  dynamic  pressure  q can  be  expressed  by 
the  equation 


q 


-at 

e 


Approximations  of  the  decay  parameter  a or  the  time  constant 

t • 1/a 


(A3) 


(A4) 


are  obtained  as  follows: 


Figure  3-21  of  Reference  11  shows  the  variation  of  dimensionless 
dynamic  pressure  q/q  with  dimensionless  time  t/t^.  The  time  t^ 
represents  the  duration  of  the  positive  phase  of  the  dynamic  pressure. 
Figure  3-21  displays  the  variation  of  q/qo  for  several  peak  overpressures 


S9MGMDUO  Pifll 


Pao  (or  peak  dynamic  pressures  qQ)  for  a 1 MI  yield  weapon.  Values  of 
t/t  read  from  the  curves  at  q/q  *0.1  are  listed  in  Table  A1  in 
column  3.  (The  subscript  is  appended  to  denote  the  time  when  q/qQ  ■ 0.1.) 
The  corresponding  psQ  and  qQ  values  for  these  curves  are  listed  in 
columns  1 and  2.  Values  of  t are  shown  in  column  4 which  have  been 
read  from  Figure  3-19  in  Reference  11,  corresponding  to  the  value  of 
p . Setting  q • 0.1  q in  Equation  (A3)  and  solving  for  a gives 

SO  o 

„ . ism 

or 

1 " in (10) 

Values  of  t listed  in  column  (5)  in  Table  A1  are  obtained  from 
Equation  (A5b) . 

Figure  A1  shows  the  variation  of  t with  qQ.  The  data  indicate 
the  variation  can  be  approximated  by  a straight  line  on  a log-log  plot. 
This  implies  qQ  can  be  approximated  by  the  equation 

inqQ  - in  (A)  + B in  t (A6) 

The  straight  line  approximation  in  Figure  A1  corresponds  to 
A ■ 0.25  and  B « -2.72,  so  that  (A6)  becomes 

inqo  --1.39  - 2.72  iu  t 


(A5a) 


(A5b) 


or 


/1.39  + inq  \ 

"\  2772  / 

x ■ e ' 


(A7) 


Equation  (A7)  represents  the  time  constant  for  a 1 MT  weapon.  The  time 
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APPENDIX  B ELASTO-PLASTlC  ANALYSIS  MODEL 

Solutions  of  the  single-degree-of-freedom  elnsto-plastic  system 
considered  are  derived  in  this  section  for  the  idealized  system 
illustrated  in  Figure  Bl. 

B.l  Elastic  Ramie 

The  governing  differential  equation  for  the  first  elastic  range 
Oiyiye  la 


-at 


my  + ky  - PQe 


y(0)  - HO)  - 0 


(Bl) 


and  the  solutions  ure 


y “ 2 2” 

m(ot  + w.  j 


[4 


— sinui^t  - cosw^t  + e 


-at 


] 


m(a- 


— 7-  lacosai^  + w.sina)^  - ae’*atl 

+ »i>  L 1 1 1 J 


(B2) 


where 


‘1  • V? 


(B3) 


If  the  disturbance  force  is  not  large  enough  to  cause  the 
displacement  to  exceed  the  elastic  limit  yg,  (B2)  are  valid  for  all 
time  in  which  case  the  displacement  increases,  reaches  a maximum  and 
then  decreases.  The  time  at  maximum  Lb  denoted  by  t2  and  is  obtained 
by  setting  y(t2)  ■ 0 in  (B2)  obtaining 


■ COSU)^  t2 


W1 

+ — 8inu).,t0 
a 12 


(B4) 
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The  corresponding  maximum  value  is  obtained  from  (B2) 


'max  , 2 
m(ot 


io U 

+“i2>  h 


-at. 


ainui^t^  - cosw^t^  + e 


which  can  be  simplified  by  making  use  of  (B4)  to  obtain 


/ ■ sinu).  t_  (B5) 

max  maw,  12 


Suppose  the  disturbance  force  is  large  enough  to  cause  the 

displacement  to  exceed  y . Then  equations  (B2)  are  valid  until  time 

t^  when  y - ye>  This  time  is  obtained  from  (B?)  by  substituting 

t » t,  and  y * y to  obtain 
1 e 

m(a2  + o)12)ye  -otj  a 

r e • — sinu),  t,  - cosw.t,  (B6) 

P w,  1 1 11 

o 1 


The  corresponding  velocity  at  time  t^  is 

P r 


-at 


y(t,) 


1 ?'T  , 2 . 2,  ' 

in(a  + ) L 


lacosw^t^  + w^sinw^t^  - ae 


-at, 


This  can  be  simplified  by  making  use  of  (B6)  to  eliminate  e 
obtain 


/,  ■ — - sinu,  t,  - ay 
1 mu^  11  Je 


■] 


to 


(B7) 


B,2  Plastic  Range 


For  those  cases  where  Pq  is  large  enough  for  the  displacement  to 
exceed  y , the  governing  equations  are  (y  > yg) 


mV  + k^y  ■ P e 


-at 


(B8) 
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The  general  solution  of  (B8)  is 


(B9) 


P e 

o 


-at 


y • A slnujt  + B cosui-t  + — —jj jT" 

m(a  + u>2  ) 


(BIO) 


where 


•V? 


(Bll) 


The  constants  A and  B obtained  by  Imposing  conditions  (B9)  are 
found  to  be 


A - C1sinw2t1  + — cosui2t1 


(B12) 


where 


B ■ C,  cosuut,  - — sinui,t. 
I Z 1 1*2  * 1 


-at- 


P e 
o 


ci  ■ ye  - 7 2 7 ~r 

m(a  + u)2  ) 


-at. 


P ae 
o 


C2  “ V1  + ,2  2. 

m(a  + u>2  ) 


(B13) 


Let  t3  represent  the  time  when  (BIO)  reaches  its  maximum  value. 
This  occurs  when  y(t3>  - 0 which  gives 
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APPENDIX  C SIDE-ON  LOADING  ANALYSIS 


C.l  Introduction 

In  this  derivation  It  Is  assumed  the  detonation  point  Is  located 
so  as  to  subject  the  side  of  the  building  to  the  blast  wave.  The  angle 
between  the  detonation  point  and  a perpendicular  to  the  rail  beams  which 
support  the  crune  is  denoted  by  ijj.  When  <t>  ■ 0,  the  detonation  point 
is  centered  on  the  side  of  the  building  and  would  not  subject  the  front 
of  the  building  to  a dynamic  pressure  wave.  When  1J1  is  different  from 
zero  both  the  aide  and  front  of  the  building  are  subjected  to  dynamic 
pressure  loading. 

Failure  i.s  considered  to  be  initiated  when  the  crane  support  beam 
deforms  sufficiently  to  ullow  the  crane  to  fall  to  the  floor  below.  The 
minimum  deformution  corresponds  to  the  crane  located  where  the  maximum 
rail  deformation  occurs,  nominally  at  the  midpoint  of  the  rail.  Failure 
ie  considered  to  be  u certuinty  if  the  deformations  are  large  enough 
to  exceed  the  plastic,  failure  limit. 

Equations  relating  weapon  yield  and  range  to  the  two  failure  limits 
are  derived  here  for  u single  degree  of  freedom  representation  of  the 
crane  rail  and  its  support  columns.  The  single  degree  of  freedom  system 
is  considered  to  be  un  clasto-plastic  system  which  may  have  a work 
hardening  behavior  after  yield  begins  in  the  rail  beam. 

Thu  lateral  displacement  of  the  crane  rail,  denoted  by  y,  corresponds 
to  the  total  deformation  that  would  be  measured  at  the  center  of  the 
rail  beam.  This  displacement  consists  of  the  sum  of  the  displacement 
of  the  top  of  the  vertical  columns  (where  the  rail  beams  are  attached) 
plus  the  displacement  at  the  center  of  the  rail  beam  relative  to  its 
end  attachment  points.  The  single  degree  of  freedom  approximation 
begins  with  a determination  of  the  linear  spring  rates  and  static  mode 
shapes.  The  spring  rate  is  determined  by  applying  a load  P at  the 
center  of  the  crane  support  beams.  One  half  of  this  load  is  reacted 
by  each  of  the  two  columns.  Equations  are  derived  for  various  end 
fixity  conditions  both  at  the  beam-column  interface  and  also  at  the 


column-ground  Interface.  Once  the  deformed  shapes  have  been  determined., 
the  generalized  mass,  spring  constant  and  aerodynamic  forces  are 
determined.  The  column-support  beam  system  and  its  single  degree  of 
freedom  approximation  are  illustrated  in  Figures  1 and  Bl.  Figure  Cl 
shows  an  idealization  of  the  elastic  members  when  subjected  to  the  load 
p.  A list  of  the  parameters  used  in  this  analysis  is  given  in  Table  C-l. 

C.2  Beam  Deformed  Shape  and  Spring  Constant 

The  moment  at  any  point  x2  (in  the  range  0 < < (£2/2)  due 

to  the  load  P is 


M ■ M - - x 
n o 2 2 


(Cl) 


The  root  moment  Mq  corresponds  to  the  moment  reacted  by  the  column. 
If  the  beam  is  simply  supported,  the  moment  Mq  is  zero.  By  assuming 
the  beam  is  uniform  one  can  write 


d \ P 

EI2  dx2  " Mo  “ 2 x2 


(C2) 


and  then  by  assuming 


dy2 

y2<*2  - 0)  - 0 and  ^ 


(-2  * t)  ' 


(C3) 


Equation  (C2)  integrates  to  give 


E*2y2 


■{4 


V2)  2 ^ 

2 / x2  2 X2  ” 12 


(C4) 


The  support  moment  Mq  can  be  eliminated  by  introducing  a fixity  factor 
F2  defined  as 


Mq  (corresponding  to  actual  system) 
E2  " M (corresponding  to  a clamped  beam) 


(C5) 
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Table  Cl  Parameters  for  Side-on  Loading  Case 


PARAMETER 

MEDIAN 

neutral  axis  to  outer  fibers  of  beam  (in) 

6.0 

column  drag  coefficient  (dimensionless) 

1.3 

beam  drag  coefficient  (dimensionless) 

1.3 

column  frontal  depth  (in) 

14.8 

beam  frontal  depth  (in) 

36 

Young's  modulus  for  column  (psi) 

3x10  7 

Young's  modulus  for  beam  (psi) 

3xl07 

column  end  fixity  factor  (dimensionless) 

0.5 

beam  end  fixity  factor  (dimensionless) 

0.5 

acceleration  of  gravity  (in/sec  ) 

386 

column  area  moment  of  inertia  (in  ) 

1700 

beam  area  moment  of  inertia  (in  ) 

300 

column  length  (ft) 

33 

beam  length  (ft) 

40 

lumped  roof  mass  (lb  sec  / in) 

25 

ambient  pressure  (psi) 

14.7 

range  to  ground  zero  (kft) 

5 

weapon  yield  (KT) 

20 

column  weight/unit  length  (lb/ft) 

140 

beam  weight/unit  length  (lb/ft) 

170 

minimum  displacement  of  beam  for  crane 
drop  (in) 

6 

work  hardening  parameter  (dimensionless) 

0.25 

column  area  factor  (dimensionless) 

2 

beam  area  factor  (dimensionless) 

2 

ductility  parameter  (dimensionless) 

6 

beam  yield  stress  (ksl) 

36 
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and  the  governing  differential  equation  is 


El 


d y:  Pi1  PXl 

1 , 2 " ~2  2~ 
dx  , 


(C13) 


Solutions  of  (C13)  subject  to  the  conditions 


V*!  " 0)  ‘ °>  (X1  ■ °)  * 81 


(CIA) 


are 


P £ p 

EIlyl  “ tEIl0l)xl  + ~ X1  " ~12 


dv  PJ,  Px.- 

EI1  d ^ " EI191  + — X!  - IT 


(C15) 


By  letting  <5^  represent  the  deflection  at  the  top  of  the  column 
(X1  ■ i^)  (C15)  gives 


6 


1 


(C16) 


Note  that  if  the  base  was  pinned,  the  system  would  be  unstable  so  one 
cannot  introduce  a fixity  factor  based  on  root  moment  that  can  be 
considered  to  be  zero  in  the  limit.  The  column  fixity  factor  is 

defined  as 


P 


1 


6 (due  to  rotation  of  base) 
(due  to  bending  of  column) 


(C17) 


With  this  definition  * 1 corresponds  to  a clamped  base,  similar 
to  ?2  ■ 1 corresponding  to  a clamped  support  for  the  beam.  However, 
P^  ■ 0 corresponds  to  a root  support  in  which  the  top  of  the  column 
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With  this  definition  F0  - 0 corresponds  to 
and  F2  ■ 1 corresponds  to  a clamped  beam. 


a simply  supported  beam. 
For  a clamped  beam 


<C6) 


Introducing  (C6)  into  (C5)  and  solving  for  Mq  gives 


M 


o 


F0  M 
2 oc 


( C 7 ) 


which  introduced  into  (C4)  gives 


El 


2y2 


(1  - f2) 


+ 


(C8) 


By  letting  6^  represent  the  deflection  at  the  center  (where  the  load 
P is  applied)  relative  to  the  top  of  the  column,  (C8)  gives 


P^3  A - 3F2\ 

°2  “ EI2  \ 192  / 

The  corresponding  beam  spring  constant  equal  to  P/62  is 

192a2 

K “ 4-  3F2 


(C9) 


(CIO) 


where 


beam 


C.3  Column  Deformed  Shape  and  Spring  Constant 


(Cll) 


The  moment  at  any  point 
to  the  end  load  P/2  is 


of  the  column  shown  in  Figure  Cl  due 


(C12) 
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has  equal  displacements  contributions  due  to  column  bending  and  column 
rotation.  By  making  use  of  (C16)  Equation  (C17)  becomes 


Mi 

Fi  • 1 r1 — 

P^' /bEl1 


Ml  “ (1  - V 6eT 


(C18) 


which  substituted  Into  (C16)  gives 


- 1 I 


S!  ■ <2  ‘ Fl>  6E t 


(C19) 


Hence,  the  column  spring  constant  k • (P/2)/<5.  becomes 

C 1 


k - t; sr 

c i — f. 


(C20) 


where 


-•W 


column 


By  making  use  of  (C19)  Equation  (C15)  gives  the  deflection  at  any 
point  to  be 

b-v^ft)  -ft)] 

C.4  Effective  Spring  Constant 

The  effective  spring  constant  of  the  single  degree  of  freedom 
approximation  is  composed  of  a parallel  combination  of  two  column 


(C21) 


(C22) 


- ■ ‘Af,  if 


sSfrtebkil/ffr’ 


springs  In  series  with  the  beam  spring.  Hence,  the  effective  spring 


constant  k Is 


<2kc)(V 

2kc+kb 


(C23) 


where  and  k^  are  given  in  (CIO)  and  (C20) . 

C.5  Relation  Between  Beam  Yield  Stress  and  Deformation 

The  deflection  y that  produces  a bending  stress  equal  to  the 
yield  stress  in  the  outermost  fibers  of  the  beam  is  denoted  by  yg. 

The  corresponding  relative  deformation  in  the  beam  is  denoted  by 
The  relation  between  the  yield  stress  a and  the  deformation  of 
single  degree  of  freedom  system  yp  is  derived  in  this  section.  The 
total  deformation  of  the  single  degree  of  freedom  system  y is  related 
to  column  and  beam  deformations  by  the  equation 


y - «!  + «2 

Since  the  same  force  ex.'.sts  in  the  equivalent  system,  the  beam  and  in 
the  columns,  one  can  write 

ky  “ ?Vl  " kb62 

From  (C23) , (C2A)  and  (C23)  one  obtains  the  relative  deformation  of 
the  beam  center  6^  and  the  deformation  of  the  top  of  the  column 
6^  for  a given  total  deformation  y to  be 


(C2A) 


(C25) 


5i'W 


(C26) 


The  maximum  bending  moment  occurs  at  the  center  of  the  beam  and  according 
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to  Equation  (C8)  is  (sign  was  changed  to  make  moment  positive) 


M 

max 


- El, 


'l2^1 


(2  - F2) 


(C27) 


The  maxinum  moment  can  be  expressed  in  terms  of  6^  by  solving 
(C9)  for  F and  introducing  the  result  in  (C27)  to  obtain 


M 


max 


24EI262 


(C28) 


and  the  corresponding  maximum  bending  stress  °max 
located  a distance  C from  the  neutral  axis  is 


in  the  outer  fibers 


o 

max 


(C29) 


At  a deflection  which  produces  the  yield  stress  one  writes  62 

and  a ■»  a and  (C29)  gives 
max  y 6 


5 


2e 


(C30) 


By  making  use  of  (C26)  the  deflection  yg  can  also  be  expressed  as 


o 

y 


24EC  (2  ~ F2  \ /k_\ 

122  V“3F27\V 


(C31) 


C,6  Generalized  Mass 

The  generalized  mass  m of  the  equivalent  single  degree  of  freedom 
system  is  obtained  from  the  equation 


2, 

y i dm. 


(y2  + Sj^)2  dm2 


(C32) 
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which  is  based  on  an  energy  equivalence  criteria.  The  factor  2 in 
the  first  integral  accounts  for  2 columns,  while  the  factor  2 on  the 
second  integral  is  introduced  because  the  integration  extends  over  one 
half  the  beam. 


In  Equation  (C32)  y represents  the  displacement  of  the  equivalent 
single  degree  of  freedom  system,  y^  represents  the  shape  of  the  deformed 
column  under  the  end  load  P/2,  y,.  represents  the  deformed  shape  of 
the  beam  under  the  central  load  P,  and  <5^  represents  the  displacement 
of  the  beam  supports.  Upon  replacing  y by  the  expressions  (C24) 

(and  then  introducing  (C19)  and  (C9)),  and  substituting  (C8)  for  y 2 
and  (C22)  for  y^  and  finally  performing  the  integrations,  the  results 
can  be  put  in  the  form 

m - K1m1  + <K2  + AK2)m2  (C33) 


where 


* .512 
K1  105 


[470  - 511^  4-  U0Fi2](^) 

[4  ' 3F2  + 32  <2  - V 

L_ 


272  - 427F2  + 168F2, 


4 - 3F2  + 32(2  - F 1 ) ^ 


K2  « 32 


" Fl>  ° j 

[4  ' 3F2  + 32(2  ' V 


C.7  Generalized  Force 

The  generalized  force  of  the  equivalent  single  degree  of  freedom 
system  is  obtained  from  the  relation 


(C34) 
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"1  2/2 
Py  - 2^  y1dP1  +2  j Cy2  + 6J 


(C35) 


which  Is  based  on  a work  equivalence  principle.  The  differential  work 
terms  arriving  from  the  dynamic  pressure  loading  are  expressed  as 

dPl  " ^VlSl5  dxl  ( 


(C36) 

(C36) 


dP2  - q(n2d2CD2)  dx2 


where  q is  the  dynamic  pressure  loading,  CQ  is  the  drag  coefficient, 
d represents  the  width  of  member  exposed  to  the  blast  pressure.  The 
n factors  are  introduced  to  account  for  a possible  increase  in  effective 
frontal  area  resulting  from  some  effect  such  as  the  outer  skin  of  the 
building  wrapping  around  the  member  and  Increasing  the  area  exposed 
to  the  air  blast.  Upon  introducing  the  expressions  for  y,  y^  and 
y2  and  performing  the  integrations,  as  was  done  for  the  generalized 
mass,  the  generalized  force  P becomes 


P - pi  + p2 


(C37) 


where 


1 2 


CU1(7  - 4F1>  ^ 


32(2  - F. ) + 4 - 3F, 

1 Ot^  l 


(C38) 
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32(2  - Fx)  ~ + 4 - 3F2 
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Figure  Cl:  Idealized  System  Elements. 

(a)  Elements  Subjected  to  Pressure  Loading 

(b)  Elements  Subjected  to  Centerload  P 


APPENDIX  D.  END-ON  LOADING  ANALYSIS 


D.l  Introduction 

In  this  derivation  it  is  assumed  the  detonation  point  is  located  such  that 
it  subjects  the  end  of  the  building  to  a blast  wave  pressure  loading.  The  angle 
between  the  detonation  point  and  a perpendicular  to  the  crane  is  denoted  by  ip. 

When  ip  * 0,  the  detonation  point  is  directly  in  front  of  the  end  of  the  building 
and  the  blast  wave  impinges  on  the  end  of  the  building  without  impinging  on  the 
side  of  the  building.  When  \p  is  different  from  zero  a component  of  the  blase 
wave  impinges  on  the  sides  of  the  building. 

From  Figure  D-l,  it  can  be  seen  that  even  if  ip  - 0 and  the  blast  wave  subjects 
the  crane  to  a uniform  pressure  loading  along  its  length,  a yaw  moment  would  exist 
about  the  mass  center  because  of  an  unsymmetrical  frontal  area  exposed  to  the 
pressure  wave.  The  rotation  produced  by  the  yaw  moment  is  resisted  by  contact 
forces  between  the  wheels  and  rail  beam.  If  the  disturbance  moment  is  large  enough 
compared  to  the  resistive  moment,  the  crane  will  rotate  and  fail  the  the  floor 
below.  Equations  governing  this  type  of  failure  are  derived  in  this  appendix. 

A list  of  the  parameters  used  in  the  derivation  are  given  in  Table  D-l. 

D.2  Single  Degree  of  Freedom  Approximation 

In  general  the  impinging  blast  wave  produces  both  a translation 
and  rotation  of  the  crane.  Translational  motion  of  even  several  feet 
may  not  necessarily  be  critical  with  regard  to  the  crane  leaving  its 
track  and  falling  to  the  floor.  On  the  other  hand,  only  a few  degrees 
of  rotation  would  be  sufficient  to  allow  the  crane  to  fall.  In  this 
section  translational  motion  is  neglected  and  only  the  yaw  rotational 
degree  of  freedom  is  considered.  By  letting  0 represent  the  yaw  angle, 
the  differential  equations  governing  the  motion  for  the  system  illustrated 

in  Figure  Dl  is 

IQB  - Mq  - £(a  - b)  F2  + bF1J  sgn  Mq 

where  I is  the  yaw  mass  moment  of  inertia  about  the  mass  center 
o 

and  F^  and  F2  are  the  rail  reaction  forces. 
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Table  D-l.  Parameters  for  End-on  Loading  Case 


PARAMETER 

MEDIAN 

crane  length  (ft) 

78 

distance  from  left  rail  to  crane  center  of  mass  (ft) 

39 

minimum  inward  motion  of  crane  before  failure  (ft) 

0.5 

maximum  inward  motion  of  crane  before  failure  (ft) 

1.0 

drag  coefficient  (dimensionless) 

1.3 

frictional  resistance  at  left  rail  (lb) 

1000 

frictional  resistance  at  right  rail  (lb) 

1000 

moment  of  inertia  about  center  of  mass  (ft  lb  sec  ) 

25000 

depth  of  crane  (ft) 

0.8 

depth  of  operator  housing  (ft) 

2.5 

distance  from  center  of  mass  to  housing  (ft) 

10 

length  of  housing  (ft) 

4 

range  to  ground  zero  (kft) 

5 

crane  weight  (lb) 

2000 

weapon  yield  (KT) 

20 

area  factor  (dimensionless) 

2 

The  blast  pressure  disturbance  moment  Mq  is 


where 


Mq  ■ q costji  / (nhCg)  udu 


M - q e ut  nCLl^cost/' 


u h iu 


5 hxa(a  - 2b)  + \ + 2^j 


In  (02)  the  change  In  the  nonnut  pressure  due  to  a rotation  is  neglected 
because  one  only  nueds  to  consider  small  changes  in  0 from  its  original 
zero  value.  Here  again  the  n factor  i«  included  to  account  for  a 
possible  increase  in  effuctive  area  upon  whi  *h  the  blast  pressure  acts. 

Equation  (01)  can  be  written  In  tlu*  form 


*• 

0 - Ac  - B 


where 


M q nil  cohiJi 
o _ o U 


M..  (a  - b)  M.  + bl'\ 

0 ' ~ * i— 1 u8*^M0) 

o o 


The  solution  of  (0A)  for  zero  itwtiul  conditions  is 
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The  time  when  6 reaches  a maximum  is  denoted  by  and  is  obtained 

by  setting  t • t^  and  0 ■ 0 in  (D6)  to  obtain 


sf1-**1)'*!- 


where 


*i  ‘ 


One  observes  that  a maximum  will  only  occur  when  A > B,  i.e., 
the  disturbance  moment  must  be  larger  than  the  resistive  moment.  When 
A > B the  maximum  yaw  angle  is 


J (“4  + e"atl  - l)'  5 “h' 


which  cun  be  written  more  simply  by  making  use  of  (D7)  and  (08)  to 
obtain 


U • 3 Two  Degree  of  Freedom  Approximation 

The  yaw  motion  0 considered  in  the  preceding  section  is  generalized 
here  to  account  for  crune  motion  along  the  rail  bourns.  Let  y^  und 
yj  represent  the  dovntrack  motion  of  the  left  and  right  ends  of  the 
crane.  One  can  describe  the  crane  motion  in  terms  of  uny  two  of  the 
three  coordinates  y^,  y^  or  0.  In  this  derivation  motion  is  described 
in  terms  of  y^  und  0. 


The  kinetic  energy  of  the  system  T is 


T " 2 m£j  ^ ■*  J + mb^jft  cosO 


(DIO) 
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where  I Is  the  yaw  mass  moment  of  inertia  about  the  left  end.  The 
dynamic  pressure  force  is  taken  to  be 


q ■ q^-01*  co*i<i>  + 0) 


and  the  track  reaction  forces  are 

Fx  - s&nCy^ 

F2  - f2  ®8n(*2) 

where 


y2  * + a tan  0 


(Dll) 


(D12) 


2 

y2  ■ y^  + a 6 sec  0 


(D13) 


The  generalized  forces  associated  with  the  disturbance  and  reaction 
forces  are 


s ; 

'? 


I 

: 


■ nbS  aln0  - coa(^i  + 6) 

G2  ■ -Fj  aaec^Q  + nCj^q^-0*  cos (tj»  + 0) 


The  solution  of  Equations  (D16)  for  and  6 gives 


Equations  (D18)  are  to  be  solved  subject  to  the  Initial  conditions 


y,(0)  - y. (0)  - 0(0)  ■ 6(0)  - 0 
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Con . 
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line . 

Con, 

for 

R, 

Schock 

ATTN: 

Dm  , 

Con. 

for 

0. 

Kahn 
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Scientific 

Laboi 

eatery 
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Doc, 

Con, 

for 

A, 

Davis 

AT  1 N : 

Doc . 

Con. 

for 

T, 

Howler 

ATTN: 

Doc. 

Con, 

for 

«. 

spt  1 linait 

ATTN: 

Doc. 

Con. 

for 

Reports  Library 

Oak  Rlduo  National  Laboratory 

Union  Carbide  Corporation-Nuclear  Division 

X-1U  Laboratory  Records  Department 

ATTN:  Doc.  Con,  fur  Civil  Uef.  Res.  t’ruj . 
ATTN:  Doc.  Con,  for  Technical  Library 

Sandia  Laboratories 
Livermore  Laboratory 

ATTN:  Doc.  Con,  for  Library  A Sucurity 
Classification  Division 
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a tori 
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Doc, 
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Doc . 

Con, 

for 
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Vortman 

ATTNl 

Hoc, 

Cun . 

for 

W. 

Rohorty 
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Hill 
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Central  Intelligence  Agency 

ATTN:  RD/SI , Rm.  5G48,  Hq.  Bldg,  for  NED/OS  I - 
5G4B.  Hqs . 

Department  of  the  Interior 

Bureau  of  Mines 

ATTN:  Technical  Library 

Department  of  the  Interior 

U.S.  Geological  Survey 
ATTN:  D.  Roddy 

NASA 

Ames  Research  Conter 
ATTN:  R.  Jackson 

Office  of  Nuclear  Reactor  Regulation 

Nuclear  Regulatory  Commission 
ATTN:  L,  Shao 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp. 

ATT N:  L.  Seizor 
ATTN:  I’.  Matliur 

Z cy  ATTN:  Technical  Information  Services 

Aybabian  Associates 
ATTN:  C.  Hague 
ATTN:  M.  Aybabian 

Analytic  Services,  Inc,  (ANSLR) 

AT1N:  C.  Ilossolbjcher 

Applied  Theory,  Inc, 
cy  A 1 1 N : J , 1 ru  I i o 

Artec  Associates,  Inc, 

ATTN:  S.  Gill 

Avco  Research  ft  Systems  Group 
ATTN:  J,  Ataimsuff 
AI1N:  Code  All  .to 
A11N:  l),  Henderson 
AI1N:  W.  11  rod  I ng 

battel le  Memorial  Institute 
ATTN:  Technical  library 
AITII:  H.  Kllngsmlth 

UDM  Corp. 

ATTN:  Corporate  Library 
ATTN:  A,  Lavagnlno 
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ATTN:  K,  Hensley 

boll  lelephone  Laboratories 
ATTN:  J,  White 

boo  ley  Co, 

ATTN:  K.  Friddell 
ATTN:  Aerospace  Library 
ATTN:  R,  Dyrdahl 
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ATTN:  >1,  Wooster 
ATTN:  It,  Holmes 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Brown  Engineering  Company,  Inc. 

ATTN:  M.  Patel 

California  'nstltute  of  Technology 
ATTN:  T.  Ahrens 

California  Research  A Technology,  Inc. 

ATTN:  K.  Kreyenhagen 
ATTN:  Technical  Library 
ATTN:  S.  Shuster 

Cal span  Corp. 

ATTN:  Technical  Library 

Center  for  Planning  A Rsch.,  Inc. 

ATTN:  R.  Shnlder 

Civil/Nuclear  Systems  Corp. 

ATTN:  R.  Crawford 

University  of  Dayton 
Industrial  Security  Super  KL-60S 
ATTN:  H.  Swift 

University  of  Denver 
Colorado  Seminary 
Denver  Research  Institute 

ATTN:  Sec.  Officer  for  J.  Wisotskl 

EGAG  Washington  Analytical  Services  Center,  Inc. 
ATTN:  Library 
ATTN:  Director 

Electric  Power  Research  Institute 
ATTN:  0.  SI  Iter 

Electromechanical  Sys.  of  New  Mexico,  Inc. 

ATTN:  R.  Shunk 

Engineering  Decision  Analysis  Co.,  Inc. 

ATTN:  R.  Kennedy 

Franklin  Institute 

ATTN:  Z.  Zudans 

Card,  Inc. 

ATTN:  (I.  Neldhardt 

General  Dynamics  Corp. 

Pomona  Division 

ATTN:  K.  Anderson 

General  Dynamics  Corp. 

Electric  Boat  Division 
ATTN:  M.  Pakstys 

General  Electric  Co. 

Space  Division 

ATTN:  M.  Bortner 

General  Electric  Co. 

Re-Entry  A Environmental  Systems  Dlv. 

ATTN:  A.  Rost 

General  Electric  Comnany-TEMPO 
Center  for  Advanced  Studies 
ATTN:  DAS  1 AC 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

General  Research  Corp. 

Santa  Barbara  Division 

ATTN:  B.  Alexander 

Geocenters,  Inc. 

ATTN:  E.  Marram 

H-Tech  Laboratories,  Inc. 

ATTN:  B.  Hartenbaum 

Honeywell,  Inc. 

Defense  Systems  Division 
ATTN:  T.  Helvig 

11T  Research  Institute 

ATTN:  Documents  Library 
ATTIi:  R.  Welch 

University  of  Illinois 
ATTN:  A.  Ang 

Institute  for  Defense  Analyses 
ATTN:  Classified  Library 
ATTN:  Director 

J,  H.  Wiggins  Co. , Inc. 

ATTN:  J.  Collins 

Kaman  AvIDyne 

Division  of  Kaman  Sciences  Corp. 

ATTN:  N.  Hobbs 
ATTN:  E.  Crlsclone 
ATTN:  Technical  Library 
ATTN:  G.  Zartarlan 

Kaman  Sciences  Corp. 

ATTN:  P.  Ellis 
ATTN:  Library 
ATTN:  F.  Shelton 
ATTN:  D.  Sachs 

Karagozlan  and  Cato 

ATTN:  d.  Karagozlan 

Lockheed  Missiles  A Space  Co.,  Inc. 

ATTN:  TiC-Llbrar,v 

Lockheed  Missiles  and  Space  Co.,  Inc. 

ATTN:  T.  Geers 
ATTN:  B.  Almroth 

Lovelace  Foundation  for  Medical  'duc^Lloe  & Rotuarth 
ATTN:  R.  Jones 

Martin  Marietta  Corp, 

Orlando  Division 

ATTN:  G,  Fotleo 
ATTN:  A.  Cowan 

University  of  Massachusetts 
ATTN i W.  Nash 

McDonnell  Douglas  Corp, 

ATTN:  R.  Italprln 

McMillan  Science  Associates,  Inc. 

ATTN:  W.  McMIlllan 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued! 

Selene#  Applications,  Inc. 

ATTN:  S.  Oston 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Merritt  CASES,  Inc. 

ATTN:  J.  Merritt 


Meteorology  Research.  Inc. 

ATTN:  U.  Green 

Mitr*  corp. 

ATTN:  Director 

University  of  New  Mexico 

Dept,  of  Campus  Security  end  Police 
ATTN:  G.  Triendeftlidis 

University  of  New  Mexico 

Civil  Engineering  Rsch.  Facility 
ATTN:  N.  Baum 
ATTN:  D.  Calhoun 

Nathan  M.  Newmark  Consulting  Engineering  Services 
ATTN:  N.  Newmark 
ATTN:  U.  Hall 

University  of  Oklahoma 

Research  Institute 

ATTN:  J.  Thompson 

Pacifica  Technology 
ATTN:  R.  Allen 
ATTN:  R.  BJork 
ATTN:  G.  Kant 

Physics  International  Co. 

ATTN:  F.  Sauer 
ATTN:  C.  Vincent 
ATTN:  R.  Swift 
ATTN:  E.  Moore 
ATTN:  D.  Orphal 
ATTN:  L.  Benrmann 
ATTN:  Technical  Library 


Science  Applications,  Inc. 

ATTN;  R.  Hoffmann 
ATTN:  D.  Maxwell 
ATTN:  D.  Bernstein 

Scienct  Applications,  Inc. 

ATTN:  U.  Layson 
ATTN:  B.  Chambers 

Southwest  Research  Institute 
ATTN:  W.  Baker 
ATTN:  A.  denial 

SRI  international 

ATTN:  W.  Wilkinson 
ATTN:  G.  Abrahamson 

Systems,  Science  A Software,  Inc. 
ATTN:  R.  Sedgewick 
ATTN:  D.  Gnn# 

ATTN:  T.  Cherry 
ATTN:  Library 
ATTN:  T.  Riney 

Terra  Tek,  Inc. 

ATTN:  Library 
ATTN:  S.  Green 
ATTN:  A.  Jones 

Tetra  Tech,  Inc. 

ATI N : L.  Hwang 
ATTN:  Library 

Texas  AIM  University  System 
c/o  Texas  A A M Research  Foundation 
ATTN:  H.  Coyle 


University  of  Pittsburgh 
School  of  Engineering 

ATTN:  M.  Williius,  Jr. 

Prototype  Development  Associates,  inc. 
ATTN:  T.  McKinley 

RAO  Associates 

ATTN:  A.  Latter 
ATTN:  P.  Rausch 
ATTN:  C.  Knowles 
ATTN:  W.  Wright,  Jr. 

ATTN:  H.  Brod# 

ATTN : J.  Carpenter 
ATTN:  J.  Lewis 

ATTN:  Technical  Information  Center 
ATTN:  A.  Fields 
ATTN:  R.  Port 


TRW  Defense  A Space  Sys.  Group 
ATTN:  A.  Narevsky 
ATTN:  N.  Lipner 
ATTN:  D.  Jortner 
ATTN:  P.  Bhutte 
ATTN:  B.  Sussholti 
ATTN:  A.  Feldman 
ATTN:  Tech.  Info.  Center 
ATTN:  J.  Ferrell 
ATTN:  0.  Chiu 
ATTN:  A.  Soux 
'i  cy  ATTN:  P.  Dai 

TRW  Defense  A Space  Sys.  Group 
San  Bernardino  Operations 
ATTN:  0.  Hulcher 
ATTN:  F.  Pieper 
ATTN:  E.  Wong 
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RIO  Associates 

ATTN:  H.  Cooper 

Rend  Corp. 

ATTN:  A.  Leupe 
ATTN:  Library 
ATTN:  C.  Mow 

Science  Applications,  Inc. 

ATTN:  Technical  Library 
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Universal  Analytics,  Inc. 

ATTN:  E.  Field 

The  Eric  H.  Wang  Civil  Engineering  Rsch,  Fee. 
ATTN:  L.  Pickle 
ATTN:  N.  Baum 

Waidllnger  Assoc.,  Consulting  Engineers 
ATTN:  J.  McCormick 
ATTN:  M.  Baron 
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Weidlinger  Assoc. , Consulting  Engineers 
ATTN:  J.  lsenberg 

Uestinghouse  Electric  Corp. 

Hiring  Division 

ATTN:  H.  Voli 


Uestinghouse  Electric  Corp. 
ATTN:  F.  Petkevicivs 


